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(57) Abstract 

A space-time signal processing system with advantageously reduced complexity. The system may take advantage of multiple 
transmitter antenna elements and/or multiple receiver antenna elements, or multiple polarizations of a single transmitter antenna element 
and/or single receiver antenna element. The system is not restricted to wireless contexts and may exploit any channel having multiple inputs 
or multiple outputs and certain otber.characteristics. Multi-path effects in a transmission medium cause a multiplicative increase in capacity. 
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SPATIO-TEMPORAL PROCESSING FOR 
5 COMMUNICATION 



10 



BACKGROUND OF THE INVENTION 

1 5 The present invention relates to digital communication and more particularly 

to a space-time communication system. 

The ability to communicate through wireless media is made difficult by the 
inherent characteristics of how transmitted signals propagate through the environment. A 
conmiunication signal transmitted through a transmitter antenna element travels along 

20 multiple paths to the receiving antenna element. Depending on many factors including the 
signal frequency and the terrain, the paths along which the signal travels will exhibit different 
attenuation and propagation delays. This results in a communication channel which exhibits 
fading and delay spread. 

It is well known that adaptive spatial processing using multiple antenna arrays 

25 increases the communications quality of wireless systems. Adaptive array processing is 

known to improve bit error rate, data rate, or spectral efficiency in a wireless communication 
system. The prior art provides for methods involving some form of space-time signal 
processing at either the input to the channel, the output to the channel, or both. The space- 
tune processing step is typically accomplished using an equalization structure wherein the 

30 time domain equalizer tap settings for a multitude of antennas are simultaneously optimized. 
This so-called "space-time equalizaion" leads to high signal processing complexity if the delay 
spread of the equivalent digital channel is substantial 

There is prior art teaching the use of conventional antenna beams or 
polarizations to create two or more spatially isolated communication channels between a 

35 transmitter and a receiver, but only imder certain favorable conditions. The radiation pattern 
cross talk between different physical transmit and receive antenna pairs must provide 
sufficient spatial isolation to create two or more substantially independent communication 
channels. This can lead to stringent manufacturing and performance requirements on the 
physical antenna arrays as well as the receiver and transmitter electronics. In addition, when 

40 large objects in the wireless propagation channel cause multipath reflections, the spatial 
isolation provided by the prior art between any two spatial subciiannels can be severely 
degraded, thus reducing communication quality. 
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What is needed is a system for more eflFectively taking advantage of multiple 
transmitter antennas and/or multiple receiver antennas to ameliorate the deleterious effects of 
the inherent characteristics of wireless media. 

5 SUMMARY OF THE INVENTION 

The present invention provides a space-time signal processing system with 
advantageously reduced complexity. The system may take advantage of multiple transmitter 
antenna elements and/or multiple receiver antenna elements, or nmltiple polarizations of a 
single transmitter antenna element and/or single receiver antenna element. The system is not 
10 restricted to wireless contexts and may exploh any channel having multiple inputs or multiple 
outputs and certain other characteristics. In certain embodiments, multi-path effects in a 
transmission medium cause a multiplicative increase in capacity. 

One wireless embodiment operates with an efficient combination of a 
substantially orthogonalizing procedure (SOP) in conjimction with a plurality of transmitter 
1 5 antenna elements with one receiver antenna element, or a plurality of receiver antenna 

elements with one transmit antenna element, or a pluraltiy of both transmitter and receiver 
antenna elements. The SOP decomposes the time domam space-time communication channel 
that may have inter symbol interference (ISI) into a set of parallel, space-frequency, SOP bins 
wherein the ISI is substantially reduced and the signal received at a receiver in one bin of the 
20 SOP is substantially independent of the signal received in any other bin of the SOP. A major 
benefit achieved thereby is that the decomposition of the ISI-rich space time channel into 
substantially independent SOP bins makes it computationally efficient to implement various 
advantageous spatial processmg techniques embodied herein. The efficiency benefit is due to 
the fact that the total signal processing complexity required to optimize performance in ail of 
25 the SOP bins is often significantly lower than the processing complexity required to jointly 
optimize multiple time domain equalizers. 

Another benefit is that in many types of wireless channels where the rank of 
the matrix channel that exist between the transmitter and the receiver within each SOP bin is 
greater than one, the combination of an SOP with spatial processing can be used to efficiently 
30 provide multiple data communication subchannels within each SOP bin. This has the 

desirable effect of essentially multiplying the spectral data efficiency of the wireless system A 
fiirther feature is the use of spatial processing techniques within each transmitter SOP bin to 
reduce radiated interference to unintentional receivers. A still fiirther feature is the ability to 
perform spatial processing within each receiver SOP bin to reduce the deleterious effects of 
35 interference fi-om imintentional transmitters. 

One advantageous specific embodunent for the SOP is to transmit with IFFT 
basis Amotions and receive with FFT basis fimctions. This particular SOP is commonly 
referred to as discrete orthogonal frequency division multiplexing (OFDM), and each SOP bin 
is thus assodaled with a frequency bin. This embodiment enhances OFDM with the addition 
40 of efficient spatial processing teclmiques. 

Accordmg to the present invention, space-frequency processing may 
adaptively create substantially independent spatial subchannels within each SOP bin even in the 
presence of significant cross talk interference between two or more physical transmit and 
receive antenna pairs. A fiirther advantage is that the space-frequency processing can 
45 advantageously adapt to cross talk interference between the physical antenna pairs even if this 
cross-talk is frequency dependent, or time varying, or both. Thus, the present invention may 
provide two or more substantially independent communication channels even in the presence 
of severe multipath and relatively poor physical antenna radiation pattern performance. 
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A further understanding of the nature and advantages of the inventions herein may be 
realized by reference to the remaining portions of the specification and the attached drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 depicts a transmitter system according to one embodiment of the 
present invention. 

Fig. 2 depicts a particular substantial orthogonalizing procedure (SOP) useful 
5 in one embodiment of the present invention. 

Fig. 3 depicts a receiver system according to one embodiment of the present 

invention. 

Fig! 4 depicts a first communication scenario where multipath is found. 
Fig. 5 depicts a second communication scenario where multipath is found. 
1 0 Fig. 6 depicts a third communication scenario where multipath is found. 

Fig. 7 depicts a multiple-input, multiple-output (MIMO) channel with 

interference. 

Fig. 8 depicts the use of an SOP in a single-input single-output (SISO) 

channel. 

15 Fig, 9 depicts the use of an SOP in a MIMO channel according to one 

embodiment of the present invention. 

Fig. 10 depicts the operation of an SOP in the context of one embodiment of 
the present invention. 

Hg. 1 1 depicts the application of spatial processing to a particular SOP bin at 
20 the transmitter end according to one embodunent of the present invention. 

Fig. 12 depicts the application of spatial processing to a particular SOP bin at 
the receiver end accordmg to one embodiment of the present invention. 

Fig. 13 depicts the application of spatial processing to N SOP bins at the 
transmitter end according to one embodiment of the present invention. 
25 Fig. 14 depicts the application of spatial processing to N SOP bins at the 

receiver end accordmg to one embodiment of the present invention. 

Fig. 15 depicts the use of a single spatial direction at the transmitter end for 
each bin of an SOP according to one embodiment of the present invention. 

Fig. 16 depicts the use of a single spatial direction at the receiver end for each 
30 bin of an SOP according to one embodiment of the present invention. 

Fig. 17 depicts the use of one or more common spatial weighting vectors for 
all SOP bins at the transmitter end according to one embodiment of the present invention. 

Fig. 18 depicts the use of one or more common spatial weighting vectors for 
all SOP bins at the receiver end according to one embodiment of the present invention. 
35 Fig. 19 depicts the use of an encoder for each SOP bin according to one 

embodiment of the present invention. 

Fig. 20 depicts the use of an encoder for each spatial direction according to 
one embodiment of the present invention. 

Fig. 21 depicts the use of an encoder for each space/fi-equency subchannel 
40 according to one embodiment of the present invention. 

Fig. 22 depicts distribution of encoder output over all space/firequency 
subchannels according to one embodiment of the present invention. 

Fig. 23 depicts a detailed diagram of an encoderAnterieaver system accordmg 
to one embodiment of the present invention. 
45 Fig. 24 depicts a transmitter system wherein multiple space/fi*equency 

subchannels are employed without spatial orthogonalization according to one embodiment of 
the present invention. 
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Fig. 25 depicts a receiver system wherein multiple space/frequency 
subchannels are employed Avithout spatial orthogonalization according to one embodiment of 
the present invention. 

Fig. 26 depicts an exemplary technique for bit loading with a trellis coder that 
uses a one-dimensional QAM symbol constellation. 



DESCRIPTION OF SPECIFIC EMBODIMENTS 

Definitions 



1 



1 0 A "channel" refers to the mput symbol to output symbol relationship for a 

communication system. A 'Vector channel" refers to a channel with a single input and 
multiple outputs (SIMO), or multiple iiqjuts and a single output (MISO). Each h entry in th e 
v ^or channel h describes one of the complex path gains present in the channel. A "matrix 
c fiSmel" refers to a channel with multiple inputs and mult i ple outputs (MIMO). Each entry 

15 Hij in the matrix H describes the complex path gain fi-om input j to output i. A "space time 
channel" refers to the input to output relationship ot a MIMO matrix channel, or a SIMO or 
MISO vector channel, that occurs when multipath signal propagation is present so that the 
channel contains delay elements that produce inter-symbol interference ^SI) as explained 
bdow. 

20 A "spatial direction" is a one dimensional subspace within a matrix or vector 

communication channel. Spatial directions need not be orthogonal. A spatial direction is 
g typically characterized by a complex input vector and a complex output vector used to weight 

— transmitted or received signals as explained herein. 

A "sub-channel" is a combination of a bin in a substantially orthogonalizing procedure 
25 (SOP) as explained below and a spatial durection within that bin. A group of spatial 
subchannels within an SOP bin may or may not be orthogonal. 

An "orthogonal dimension" is one member in a set of substantially orthogonal spatial 
directions. 

A channel "subspace" is a characterization of the complex m-space direction occupied 
30 by one or more m-dimensional vectors. The subspace characterization can be based on the 
instantaneous or average behavior of the vectors. A subspace is often characterized by a 
vector-subspace of a covariance matrix. The covariance matrix is typically a time or 
frequency averaged outer product of a matrix or vector quantity. The covariance matrix 
ch^^erizes a collection of average channel directions and the associated average strength 
35 for each direction. 

A *two norm" metric for a vector is the sum of the squared absolute values for the 
elements of the vector. 

A "Euclidean metric" is a two norm metric. 

"Intersymbol interference" (ISI) refers to the self interference that occurs between the 
40 delayed and scaled versions of one time domain symbol and subsequent symbols received at 
the output of a delay spread communication channel. The channel delay spread is caused by 
the difference in propagation delay between the various multipath components combined with 
the time domain response of the RF and digital filter elements. 

A "substantially orthogonalizing procedure" (SOP) is a procedure that plays a part in 
45 transforming a time domain sequence into a parallel set of substantially orthogonal bins, 

wherein the signals in one bin do not substantially interfere with the signals fi'om other bins. 
Typically, the transformation fi'om a tune domain sequence to a set of substantially 
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orthogonal bins requires a transmitter SOP witli a set of input bins, and a receiver SOP with a 
set of output bins. 

"Convolutional bit mapped QAM" (CBM-QAM) is the coding system that results 
when the output of a convolutional encoder are grouped and mapped to QAM constellation 
points. 

Fading "structure" occurs when the fading behavior of one or more entries in a 
channel matrix within an SOP bin is correlated across time, or frequency, or both. This 
structure can be exploited using advantageously designed estimation filters to improve 
channel estimation accuracy given multiple frequency samples of the channel matrix entries, 
or multiple time samples, or both. 

A "maximum likelihood sequence detector" is a sequence estimator that computes the 
most likely transmitted code sequence, from a set of possible sequences, by minimizing a 
maximum likelihood cost function. 

An "antenna element" is a physical radiator used to transmit or receive radio 
frequency signals. An antenna element does not involve any electronics processing 
components. A single radiator with two polarization feeds is viewed as two antenna 
elements. 

An "antenna array" is a collection of antenna elements. 

A "burst" is a group of transmitted or received communication symbols. 

Background Material 

The disclosure herem assumes a background in digital conmiunication and linear 
algebra. The following references are incorporated herein by reference. 

Wozencraft & Jacobs, Principles of Communication Engineering (1965). 

Haykin, Adaptive Filter Theory, 2"^ Ed. (1991). 

Strang, Linear Algebra. 3"* Ed. (1988). 

Jakes, Aficrowave Mobile Conununication (1974). 

Proakis, Digital Communications (1995). 

Transmitter Overview 

Fig. 1 depicts a transmitter system in accordance with one embodiment of the present 
invention. Typically an information signal input 2 mcludes a digital bit sequence, although 
other forms of digital data or analog data are possible. In the case of digital data, the input 
data sequence is first fed into an Encoder and Interleaving apparatus 10 where the data is 
encoded into a symbol stream. The symbol stream is typically a sequence of complex 
digitized values that represent members of a finite set. Each symbol can be a one dimensional 
value, or a multidunensional value. An exemplary one dimensional symbol set is a PAM 
constellatioa Note that in this discussion, it is understood that a symbol with in-phase and 
quadrature components, is considered to be a complex one dimensional symbol, so that the 
QAM constellation is also viewed as a set of one ^ensional symbols. An example 
multidimensional symbol set is a sequential grouping of QAM constellation members. 

The purpose of the encoding process is to improve the bit error rate of the transmitted 
signal by introdudng some form of mformation redundancy into the transmitted data stream. 
Usefiil encoding techniques can involve combinations of a number of well known techniques 
such as convolutional encoding with bit mapping to symbols, trellis encoding, block coding 
such as cyclic redundancy check or Reed Solomon coding with bit mapping or Automatic 
Repeat Queing. An interieaver is often advantageous for distributing the transmitted 
information among the various subchannels available for transmission. This interleaving 
distributes the effects of channel fading and interference so that long sequences of symbols 
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with poor quality are not grouped closely together in the SOP bin sequence that is fed mto 
the receiver decoder. In many applications, it is advantageous to perform a power and bit- 
loading optimization wherein the number of bits that are mapped to a given encoder symbol, 
and the signal power assigned to that symbol, are determined based upon the measured 
5 communication quality of the space-frequency information subchannel that carries the symbol 
stream. 

After the digital data is encoded into a sequence of symbols, a Training Symbol 
Injection block 20 may be used to place a set of known training symbol values in the 
transmitter symbol stream. The purpose of the training symbols is to provide a known input 

10 within a portion of the transmitted symbol stream so that a receiver may estimate the 

conununication channel parameters. The channel estimate is used to aid in demodulation and 
decoding of the data sequence. The training symbols may be injected periodically in time, 
periodically in frequency, or both. It will be obvious to one skiUed in the art that blind 
adaptive spatial processing techniques can be utilized within each SOP bin at the receiver as 

15 an alternative to training with known symbols. In such blind detection implementations. 
Training Symbol Injection block 20 is unnecessary. 

The data plus training symbol stream is then fed into a Transmitter Space-Frequency 
Pre-Processor (TSFP) block 30. The TSFP block 30 performs two sets of advantageous 
processing steps on the symbol stream before transmission. One processing step 

20 accomplished within the TSFP is the transmitter portion of a substantially orthogonalizing 
procedure (SOP). When the transmitter portion of the SOP is combined with the receiver 
portion of the SOP, a set of parallel bins are created in such a manner that information 
> ti^4^ transmitted within one bin does not substantially interfere with information transmitted from 
another bin after the receive portion of the SOP is completed. One preferred SOP pair is the 

25 inverse fast Fourier transform (IFFT) at the transmitter combined with the FFT at the 
receiver. Another advantageous SOP pair embodiment is a bank of multiple filter and 
frequency converter pairs (multi-band SOP) with one filter bank located at the transmitter and 
one filter bank located at the receiver as depicted in Fig. 2. Several other example SOPs 
including the Mbert transform pair and generalized wavelet transform pairs v^U be obvious 

30 to one skilled in the art. The other processing step accomplished in the TSFP is spatial 
processing. The spatial processing step typically multiplies one or more symbols that are 
destined for transmission in a given SOP bin with one or more spatial vector weights. For 
convenience in the following discussion, the collection of spatial processing weights applied 
to the signals transmitted or received in a given SOP bin are sometimes referred to as a 

35 matrix. The spatial vector weights are optimized to obtain various desirable performance 
enhancements. 

Fig. 2 depicts a digital baseband filter bank at the receiver and a filter bank located at 
the transmitter. Each filter of the transmitter filter bank includes a nuxer (fi^quency 
converter) 60, a bandpass filter 70 and an interpolator 80. Each filter of the receiver filter 

40 bank includes a bandpass filter 90, a nuxer 60, and a dedmator 1 00. 

One transmitter embodiment optimizes the transmitter spatial vector weights so that 
the multiple subchaimels in a given SOP bin can be converted at the receiver into substantially 
independent received spatial subchannels wherein symbols from one subchannel do not 
substantially interfere with symbols from another subchannel. Another embodiment optimizes 

45 the transmitter spatial vector weight to improve the received power of one or more spatial 
subchannels within each SOP bin, or to improve the average power of several spatial 
subchaimels within several SOP bins. A fiuther embodiment optimizes the transmitter spatial 
vector weights within each SOP bin to simultaneously increase the power delivered to the 
desired receiver within one or more spatial subchannels while reducing interference radiated 
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to unintended receivers. A yet further embodiment spatial processes one or more symbols 
within each SOP bin by multiplying each symbol with a transmitter weight vector that is fixed 
for all SOP bins, with the weight vectors optimized to increase the time or fi-equency average 
power delivered to one or more desired receiver spatial subchanels, and possibly reduce the 

5 time or fi-equency averaged interference radiated to unintended receivers. This last 

embodiment is particularly useful in FDD systems where miiltipath fading makes it impossible 
to esthnate the forward channel fi-om reverse chaimel data, but where the average forward 
channel subspaces are substantially similar to the average reverse channel subspaces. Another 
embodiment teaches simply routing each symbol firom the encoder to one antenna element in 

10 each SOP bin without any weighting. Other usefijl embodiments are discussed herein, and 
many others usefiil combinations of spatial processing with an SOP will become obvious to 
one skilled in the art. It is understood that one or more digital filters are typically used in 
TSFP 30 to shape the transmitted RF signal spectrum. 

Once the encoder symbol sequence is processed by TSFP 30, the processed symbol 

15 sequence includes a parallel set of digital time domain signal sequences. Each of these time 
domain signal sequences is fed into one input of a Modulation and RF System block 40. 
Modulation and RF System block 40 uicludes a set of independent RF upconverter chains that 
fi-equency convert the digital baseband signal sequence up to the RF carrier frequency. This 
is accomplished using apparatus that mcludes digital to analog converters, RF mbcer 

20 apparatus, and fi*equency synthesizer apparatus. The details of these elements of the 
invention are well known and will not be discussed here: 

The final step in the transmission process is to radiate the transmitted signal using a 
l^ifi^ii Transmit Antenna Array 50. The antenna arrays can be constructed from one or more co- 
polarized radiating elements or there may be multiple polarizations. If there is multipath 

25 signal propagation present in the radio link, or if there are multiple polarizations in the 
antenna arrays, or if at least one of the antenna elements on one side of the link are in a 
disparate location from the other elements on the same side of the link, then the invention has 
the advantageous ability to create more than one subchannel within each SOP bin. It is 
understood that one physical antenna reflector with a feed that has two polarizations is 

30 considered as two antenna elements in all that follows. There are no restrictions on the 
antenna array geometry or the geometry of each radiating element. A transmitter system 
invention may adapt to provide optimized performance for any arbitrary antenna array. 

Receiver Overview 

35 Fig. 3 depicts a receiver system according to one embodiment of the present 

invention. The RF signals from each of the elements of an Antenna Array 1 10 are 
downconverted to digital baseband using a Demodulation and RF System 120. Demodulation 
and RF Systeni 120 inchides the RF signal processing apparatus to downconvert the RF 
carrier signal to a baseband IF where it is then digitized. After the digitizer, a timing and 

40 frequency synchronization apparatus is used to recover the timing of the transmitted digital 
signal sequence. Several known techniques may be used for the purpose of synchronization 
and these techniques will not be discussed herein. 

In certain embodiments of the invention, after Demodulation and RF System 120, the 
digital baseband signal is then fed into a Channel ED block 130 and a Receiver Space- 

45 Frequency Processor (RSFP) block 140. Within Channel ID block 130, the characteristics of 
the digital conmiunication channel are estimated. The estimated channel values consist of 
entries in a matrix for each SOP bin. The matrix contains complex values representing the 
magnitude of the spatial channel within the SOP bin from one transmit antenna element to one 
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receive antenna element. The matrix channel estimate for each SOP bin is provided to RSFP 
block 140 and Decoder and Deinterleaving block ISO. 

Some embodiments of the invention involve improving channel estimation 
performance by exploiting the structured nature of the frequency domain feding that exists in 
5 the matrix channels across SOP bms, exploiting the structure in time domain &ding of the 
matrix channels, or exploiting both the frequency and time domam fading structure that is 
present. By exploiting the frequency domain fading correlation, the entire set of matrix 
channels within the SOP bins may be estimated even when tr^iing information is transmitted 
in a subset of the SOP bins. This allows for simultaneous transmission of training and data 

1 0 thus reducing overhead. By exploiting the time domain correlation of the channel feding 

within each SOP bin, channel estimation accuracy is increased for a given time epoch between 
traming events. This reduces the required frequency of training symbol transmission and thus 
further reduces training overhead. It is understood that it is also possible to separately exploit 
time domain and frequency domam correlation, with the most beneficial results occurring if 

1 5 both correlation dimensions are used advantageously. It is to be understood that Chatmel ID 
block 130 is shown as a separate function even though it may share some elements with 
RSFP block 140 or Decoder and Deiriterleaver block 150. 

RSFP block 140 performs the receiver signal processing that is the dual of the two 
sets of operations performed in TSFP 30. One of the steps performed in RSFP 140 is the 

20 receiver half of the SOP. As discussed above, the receiver half of the SOP completes the 
transformation between the time domain chaimel with ISIto the substantially orthogonal set 
of bins. The second set of signal processing operations that can be performed in the RSFP is 
<J^^^^ spatial processing. In one class of embodiments, the receiver spatial processing step 

combines the output of the SOP bins using one or more vector weighted inner product steps 

25 to form one or more one-dimensional received spatial subchannels within each SOP bin. The 
receiver wdght vectors are chosen to optimize an advantageous performance measure. In 
one embodiment, wherein both the transmitter and receiver have knowledge of the channel 
state information within each SOP bin, the transmitter spatial weight vectors and the receiver 
spatial weight vectors are both chosen to optimize performance for a set of multiple 

30 substantially independent subchannels within each SOP bin. As discussed above, this 

significandy increases the spectral efiEtdency of the system. In another embodiment wherein 
the transmitter does not have channel state information, the receiver performs the spatial 
processing required to create multiple substantially independent subchannels within each 
SOP bin. In a further embodiment wherein the transmitter may or may not have channel state 

3 5 information, the receiver reduces the effects of interference radiated from unintentional 
transmitters as well as performmg the spatial processing required to create multiple 
substantially independent subchannels within each SOP bin. A yet further embodiment 
optimizes the receiver spatial vector weights within each SOP bin to simultaneously increase 
the received power and reduce the detrimental effects of interference received from 

40 unintentional transmitters. An additional embodiment involves forming one or more vector 
weights, that are fixed for aU SOP bins, where the vector weights are optimized to 
simultaneously increase the time or frequency averaged received power for one or more 
spatial subchannels, while possibly also reducing the time or frequency averaged interference 
power received from unintended transmitters. 

45 As discussed herein, certain embodiments involve simply passing the vector samples 

received in each SOP bin to Decoder and Deinterleaving block 150 without performing any 
spatial processing. It will be obvious to one skilled in the art that other combinations of 
transmitter spatial weight vector optimization techniques and receiver spatial weight vector 
optimization techniques can be constructed around the principle concept of spatial processing 
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in combination with an SOP. Other embodiments are discussed herein. One experienced in 
the art will be able to recognize additional embodiments that involve advantageous 
combination of an SOP with spatial processing at the receiver or the transmitter. It is 
understood that one or more digital filters are typically used in RSFP block 1400 to sh^pe the 
5 received RF signal spectrum. 

The outputs of RSFP block 140 are fed into Decoder and Deinterleaving (DD) block 
150. There are two broad exemplary classes of operation for the DD block 150. In the first 
exemplary broad class of embodiments, DD block 150 decodes a symbol sequence which 
was encoded and transmitted through a multitude of SOP bins with one or more substantially 

10 independent subchannels. The decoder includes the appropriate receiver counterparts for the 
combination of encoders selected for the transmitter. A preferred embodiment includes a 
deinterleaver, a trellis decoder or convolutional bit mapping decoder employing a scalar 
weighted Euclidean maximum likelihood sequence detector, followed by a Reed Solomon 
decoder, followed by an ARQ system to correct Reed Solomon codeword errors. In the 

1 5 second exemplary broad class of embodiments, DD block 1 50 decodes a sequence of 
multidimensional symbols, or groups of adjacent one dimensional symbols, with each 
multidimensional symbol or group of one dimensional symbols being received in an SOP bin. 
Typically, the symbol sequences are transmitted without weighting or with weighting that 
optimizes some measure of average signal quality. 

20 In an alternative embodiment, trellis encoded symbols are grouped and interleaved in a 

manner such that the symbols transmitted fi-om the antenna elements within a given SOP bin 
form a vector that is drawn fi"om either a multidimensional QAM encoder output symbol, or a 
J sequence of one dimensional QAM encoder output symbols that have adjacent locations in 

the pre-interleaved encoder output sequence. In this way, a maximum likelihood vector 

25 decoder may be constructed given an estimate of the ch^el matrix that is present within 
each SOP bin. The maximum likelihood decoder computes the weighted vector Euclidean 
metric given the deinterleaved received vector fi-om each SOP bin, the deinterleaved matrix 
channel estimates firom each SOP bin, and the transmitted vector symbol trellis state table. 
In another alternative embodiment, either of the aforementioned encoder 

30 embodiments will have preferable performance if the encoder polynomial and symbol 
constellation set are optimized to improve the bit error rate performance given the 
characteristics of the matrix charmel fading that occurs in each SOP bin. One particular 
metric that is well suited for a code polynomial optimization search is the product of the two 
norms of the vector diflference between the correct transmitter symbol vector and the error 

35 symbol vector. 

The output of DD block 1 50 is the estimated bit stream at the receive end of the 
radio link. 

It is to be understood that all transmitter embodiments of the present invention may be 
adapted for use with a receiver accessing the channel through a single channel output such as 
40 a single receiver antenna element. Furthermore, all receiver embodiments of the present 
invention may be adapted for use with a transmitter accessing the channel through a single 
channel input such as a single transmitter antenna element. It is understood that the channel 
is then a vector channel. Such multiple-input single-output (MISO) and smgle-input multiple 
output (SIMO) systems are within the scope of the present invention. 



45 



Space-Frequency Communication in a Multipath Channel 

Before developing the signal processing of the present invention, a physical 
description and mathematical description of wireless channels are provided. Many wireless ^ 
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communication channels are characterized by multi-path, where each path has associated 
fading and propagation delay. Multipath may be created by reflections in the physical rac&o 
path, multiple antenna polarizations, antenna elements located m disparate locations, or a 
comtsination of any of these. One scenario in which multipath is created is illustrated in Fig 4. 
5 A Base 1 52 transmits information to and receives information from a remote xinit 1 70 A or 
170B. Base 152 possesses one or more antenna elements referred to as an array 55,. 
Similarly, the Remote Units possess their own arrays 55. A transmitted signal propagates 
along multiple paths 155A-C created by reflection and scattering from physical objects in the 
terrain 160A-D. 

10 Multipath signal propagation such as that depicted in Figs. 4-6 can give rise to 

spatially selective fading, delay spread, frequency &ding, and time fading. Spatial &ding 
occurs as the various wavefronts arriving at the receiver from diflferent propagation paths 
combine with coiistructive and destructive interference at different points in space. An 
antenna array located within this spatially selective field will sample the field at various 

15 locations so that the signal strei^h at each array element is different. Delay spread occurs 
due to the differing propagation path lengths. The channel delay spread gives rise to a 
frequency selective digital communication channel at each antenna element. This frequency 
response is different for each antenna element by virtue of the fi^quency dependent spatial 
fading. Finally, if either the transmitter, or the receiver, or objects in the terrain are moving, 

20 the frequency selective spatial fading will vary as a fimction of time. The present invention is 
unique m that it is capable of eflBciently and economically adapting to the time varying space- 
frequency channel response to make advantageous use of the inherent properties of such 
^fi^iZi channels. 

For several decades, the primary focus of the prior art has been to somehow mitigate 

25 the effects of the multipath channel. This conventional approach is ill-advised since multipath 
channels give rise to a multiplicative capacity effect by virtue of the fact that the multipath 
induces a rank greater than one in the matrix channel present in each SOP bin. This provides 
opportunity to form multiple parallel subchaimels for communication within each SOP bin. 
Thus, one should utilize multipath to improve communication performance rather than 

30 attempt to mitigate its effects. A substantial advantage provided by the present invention is 
the ability to efficiently and economically exploit the inherent capacity advantages of 
multipath channels using a combination of an SOP and spatial processing or spatial coding. 
No other structures are known to efiBdentiy exploit this fimdamental advant£^e in the 
presence of substantially frequency selective multipath channels. 

35 Fig. 5 illustrates another wireless chaimel scenario in which multipath is present. A 

Base 152 with an antenna array 55 transmits information to and receives information from a 
Remote 170C with an antenna array 55. In this case, both Base and Remote antenna arrays 
55 both have elements with differing polarizations. Thus multipath signal propagation exists 
even if there are no significant reflections in the physical environment. The direct line of sight 

40 paths 155B, 155E each corresponding to one of the polarizations in the array elements, are 
suflBcient to create a matrix channel with a rank greater than one within each SOP bin, even if 
the other reflected radio paths 155A, 155C, 155D, and 155F are insignificant or nonexistent. 

It is known that such line of sight polarization (reflection free) channels can be 
decomposed mto two parallel communication channels by using high performance dual 

45 polarization antennas, one at each end of the radio link, in combination with high performance 
receiver and transmitter electronics apparatus. In this prior art, the quality of the parallel 
communication channels is limited by the degree to which the two polarization channels 
remain independent. In general, maintaining the manufacturing tolerances and installation 
alignment precision in the antennas and electronics required to achieve substantially 
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orthogonal spatial subchannels at the output of the physical receive antenna is relatively 
expensive. Slight manufacturing errors and component variations can lead to a significant 
cross-talk interference between the multiple polarizations present in the radio channel. In 
contrast, an advantageous feature of the invention is that ^e different polarizations present in 
5 the wireless channel may have an arbitrary degree of cross-talk interference, and the cross- 
talk mterference may be firequency dependent without loss of performance. In such cases, the 
invention provides an economical and efiSdent method fbr fiilly exploiting the multi- 
dunensional nature of the multiple polarization channel. It is understood that the invention 
can provide fiirther capacity advantage if the multiple polarization channel also has reflected 

10 signal paths. This additional multipath results in an additional increase in the channel matrix 
rank in each SOP bin that can be fUrther exploited to improve the capacity of the channel. 

Fig. 6 depicts another wireless commimication scenario in which multipath is present 
and can be explohed to create multiple dimensions for communication. In Fig. 6, two Bases 
152A and 152B with antenna arrays 55 communicate with Remote Units 170 A and 170B that 

1 S also possess antenna arrays 55. In this case, the composite channel is defined as the MIMO 
channel between the antennas of the two Bases 152A and 152B and the antennas of the 
Remote Units 170A and 170B. Note that this channel includes direct line of sight paths 155B 
and 153B as well as the reflected paths 155A, 155C, 153A, 153C, and 153D. By virtue of 
the spatial separation between two Bases 152 A and 152B, even if the reflected paths are 

20 insignificant or nonexistent, this channel contains multipath that can be exploited using the 
invention. In addition, the channel fi-om the antennas of the two Bases to the antennas of one 
Remote is again a matrix channel, with rank greater than one, within each SOP bin so that 
' ^^-f^^ multiple parallel dimensions for transmission may be created. In these types of applications, 
the present invention provides for the abiUty to reduce interference radiate to unintentional 

25 receivers. Furthermore, the present mvention provides the capability of reducing the 
detrimental effects of received interference fi'om unintentional transmitters. 

Thus it can be seen that multiple transmitter antenna elements or multiple antenna 
elements may be either co-located or be found at disparate locations. 
The following symbol chaimel model applies to all of the above multipath radio propagation 

30 cases illustrated by Figs. 4-6. The channel impulse response includes the effects of ihe 
propagation environment, as well as the digital pulse shaping filters used m TSFP 30 and 
RSFP 140 , the analog filters used m Modulation & RF System 40 and Demodulation & RF 
System 120. Due to the difference in propagation delay between the various multipath 
components combined with the time domain response of the RF and digital filter elements, a 

35 smgle symbol transmitted into the channel is received as a collection of delayed copies. Thus, 
delayed and scaled versions of one symbol interfere with other symbols. This self interference 
effect is termed intersymbol interfi^-ence or ISI. The delay spread parameter , denoted by v, 
is the duration in symbol periods of the significant portion of the channel impulse response. 
As the transmitted symbol rate is increased or as the physical geometries in the 

40 channel become larger, the delay spread can become so large that conventional space-time 
processing systems become highly complex. An advantage of the present invention is that the 
signal processing complexity is relatively low even whm the delay spread become s extremely 
jarge^ Thi s allo ws for the e conomical application of MIMO space-frequency proceS^g 
techniques at high digital data rates. This efiBcient use of sigmd processmg comes about 

45 because the invention allows the spagg jime channel to be t reated as a set of substantially 
independent spatials ubchannels y^thomlacrificing c hanndcapaci^ In contrast, 
conventional approaches-eitheLgttCTi pt to equaliz e the much more complex space-time 
channel or alternatively sacrifice capacity. 
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The channel is modeled as time-invariant over the time spanned by a burst of N data 
symbols, but varying from one burst to another. This block time invariant assumption 
produces a channel model that is sufBciently accurate for channels wherein the block duration 
is short compared to the channel fading, or {N + 2v)r « A^, where is the correlation 

S interval\ Other models are available wherein the channel varies continuously, but these 
models add unnecessary complication to the present discussion. It is understood that r^id 
time variation in the channel can be another motivation for choosing one of the other SOP 
alternatives in the presence of fading rates that are rapid with respect to the burst frequency. 
One skilled in the art will be aware of the pertinent issues for a given application. For 

10 example. Orthogonal Frequency Domain Multiple?dng (OFDM) is an SOP composed of an 
IFFT and cyclic prefix as the transmitter SOP, and an FFT as the receiver SOP. With OFDM, 
one pertinent issue is frequency domain inter-carrier interference (ICI), which can occur in 
OFDM systems with extremely rapid fiiding. Such pertinent issues shape the appropriate 
choice of channel models for various SOP basis fimctions. 

1 5 With this background discussion, it can be verified by one skilled in the art that the 

relationship between the transmitted burst of baseband symbols and the received burst of 
baseband samples may be adequately expressed as the space-time equation. 



x(Ar)=:G(A)z(*) + I(*), 
20 where the index k represents bursts. The composite channel output for one burst of data, 

x(k) , is written with all tune samples appear in sequence for every receive antenna 1 to A/^ , 

x(fc)=:|x,(l) x,(N^v-l) - x^(l) x^(iV+v/-l)f. 
Likewise, the input symbol vector is written, 

z(*) = |z,(l) - z,(N) ... r^^O) - ^Mri^f- 
25 The quantity I(k), defined the same as x(k) and z(k), represents both noise and interference. 
The MIMO channel matrix, G(k) , is composed of sin^e-input single-output (SISO) sub- 
blocks. 



Furthermore, each of the SISO sub-blocks, G,^ (k) , is a Toeplitz matrix describing the input- 

30 output relationship between the transmitted symbol burst and the received symbol burst for 
anterma pair (J J). This MIMO space-time chamiel is illustrated by Fig. 7, which shows SISO 
sub-blocks 180A-D and the addition of interfierence for each receiver sample. 
Space-Frequency Processing 

Embodiments of the present invention uses space-fixquency processing at either the 

3 5 transmitter or receiver, or both, to create effective communication systems in wireless 
channels. Generally, the processing substantially eliminates the ISI caused by the channel 
correlation across space (antenna correlation) and time (delay spread). This processing 
greatiy shnplifies the design of the remaining functions that comprise a complete 



^ The correlation interval here is defined as the time period required for the 
fading parameter time-autocorrelation function to decrease to some fraction of the zero-shift 
value. 
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communication system, including coding and modulation. Furthermore, the processing 
approach is based upon a capacity-achieving structure for the MIMO wireless commimication 
channel. Space-frequency processing is composed of one or more of the following: an SOP, 
a transmit spatial processor, and a receive spatial processor. 

5 

SubstantiaOy Orthogonalizing Procedure 

The use of an SOP in a SISO channel is considered first in order to illustrate the 
invention's ability to eluninate ISI across space and time. The SOP is composed of signal 
processing operations implemented at both the transmit and receive sides of the channel. This 

10 is illustrated in Fig. 8 where a Transmitter SOP processor 190 and a Receiver SOP processor 
200 jointly perform a complete SOP. The SOP ensures that the N input symbols, in bin 1 
through bin N, are transmitted through the channel in such a way that each output symbol is 
substantially influenced by only the input symbol of the same frequency bin. For example, the 
input symbol in bin 1 is the only symbol to have substantial influence on the output symbol 

15 value in bin 1. 

This concept generalizes to the MIMO system as shown in Fig. 9. For the MIMO 
system, each transmitter antenna 51 is preceded by one of Mr identical Transmitter SOP 
processors. Likewise, each receiver antenna 1 1 1 precedes one of Mr identical Receiver SOP 
processors. Hence, the processing path for any transmitter-receiver antenna pair contains a 

20 jomtly performed complete SOP. In other words, there exist MrMt SISO SOPs in the 
MIMO system. By exploiting the property of superposition, this collection of SISO SOPs 
comprise a MIMO SOP where any two symbols communicated in diflFerent bins exhibit 
' ^^^^ substantially reduced crosstalk, irrelevant of the antennas by which the symbols were 

transmitted and received. Therefore, the SOP establishes N substantially independent MIMO 

25 spatial channels. 

Many different SOP implementations exist, including an IFFT-FFT pair, a bank of 
multiple narrow-band filters, and generalized wavelet transform pairs. One advantageous 
example of an SOP is the use of a frequency transform combined with a burst cyclic prefix 
application procedure processor 207, as shown in Fig. 10. There is also a cyclic prefix 

30 removal procedure processor 206 at the receive end. When the fi^uency transform is an 
IFFT-FFT pair 205 and 208 as shown in Fig. 10, this particular SOP is commonly referred to 
as discrete orthogonal frequency division multiplexing (OFDM). Hence, this embodhnent of 
the invention combines the OFDM SOP with multiple input antennas, or multiple output 
antennas, or both multiple input and multiple output antennas. The present embodiment is 

3 5 thus termed matrix-OFDM (MOFDM). 

The analysis presented for MOFDM will have a substantially similar form as other 
choices for the SOP. These alternative embodiments have certain advantages and drawbacks 
as compared to the OFDM SOP. For example, the multi-band SOP does not completely 
eliminate ISI, but it is more robust to certain types of narrow-band interfering signals because 

40 the interference can be more confined within a given SOP bin as compared to the OFDM 
SOP. The ISI that can be present in the multi-band SOP could make it advantageous to use 
pre-equaiization or post-equalization structures in conjunction with the spatial processing 
within a given SOP bin. While this complicates the spatial processing, the complexity 
drawback may be outweighed by other requirements such as robustness to interference or the 

45 need to separate the SOP bins by relatively large frequency separation. Only the OFDM SOP 
will be analyzed m detail herein and it will be understood that one may exploit the other SOP 
choices as needs dictate. 

As depicted in Fig. 10, the exemplary SOP operates in the foUowing fashion. The 
symbols from the transmit spatial pre-processor, z. (n) , considered to be in the frequency 
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domain, are organized into Mj vectors of N complex symbols. Each of these vectors is then 
converted to the time-domain using aniV^point inverse-fast-Fourier-transform (BFFT) 
procedure 205. Each of the parallel My. time-domain sequences has a cyclic prefix added to 
the beginning, so that the last v elements in the IFFT output sequence form a pre-amble to 
5 the JV^element IFFT output. The cyclic prefix operation is given by; 

[z(\) z{N)Y 1-4 [r(J\r-i/ + l) z{N) z(l) - z{N)Y 
The application of the cyclic prefix is performed by cyclic prefix application procedure 
processor 207. For each antenna, the ( iV^ + v )-length sequences are passed to the RF 
transmit chain for D/A conversion and modulation. 

10 Likev^se, each RF receiver chain produces a sampled sequence of length N+v, 

Cyclic prefix removal procedure processors 206 remove the cyclic prefix fi'om each sequence 
by discarding the first v data symbols, resulting in Mj^ vectors of complex symbols. Each 
of these sequences is then processed with an N-point fast-Fourier-transform (EFT). 
These symbols are then passed to the receiver spatial processor. 

15 The eflFect of the SOP is to substantially remove the ISI between any two symbols 

assigned to different bins, for any pair of transmit and receive antennas. Therefore, for each 
IFFT-FFT bin «, the received signal values for each antenna, x(w) , are related to the 
transmitted firequency-domain symbols, z(ff), through the expression, 

x(w) = H(/i)z(n)+I(/i) Vw (10) 

20 where x(n) is a complex -element vector at SOP bin «, z(n) is a complex Mj -element 
^ symbol vector at bin /i, and I(n) is the interference and noise at all receive antennas for bin n. 
Note that a time index is not included in the above equation since it is assumed that cbatmel is 
time-invariant over the length of a burst. The spatial sub-channels, H(w) , are M^^ by Mj 
element matrices that describe the spatial correlation remaining m the wu*eless channel after 

25 the SOP. For the MIMO case, each SOP bin may be characterized by a matrix of complex 
values, vwth each value representing the path gain fi-om a given transmit antenna element to a 
given receive antenna element in that particular SOP bin. 

To understand the result given by Equation(lO), it is instructive to show how the SOP 
pre-processor and post-processor acts upon the time-domain channel. The MIMO time- 

30 domain channel, G<k), contains Mj^Mj Toeplitz matrices that describe the time-domain 
iiqjut-output behavior of each antenna pair (see eq. 1). This channel formulation is depicted 
in Fig. 10. It is well known that by adding a cyclic prefix at the transmitter and subsequently 
removing the prefix at the receiver, a Toeplitz intput-output matrix is transfomed into a 
circulant input-output matrix (the /th row is equal to the yth row cyclicly shifted by i-j 

35 elements). Therefore, the each G, ^ m Fig. 10 is transformed into a circulant G, ^ . The 

MIMO circulant matrices are delimited in Fig. 10 by G . 

This particular class of SOP exploits the fact that any circulant matrix can be 
diagonalized by a predetermined matrix operator. One such operator is a matrix of the FFT 

basis vectors. That is, for any caculant matrix G , 
40 D = YGY'' 

where Dis some diagonal matrix and the scalar elements of Y are. 

Applying Mj IFFT "operations at the transmitter and FFT operations at the receiver is 
described mathematically by a pre-multiplication of a NM^ x NM^^ block diagonal FFT matrix 



wo 98/09385 



PCTAJS97/15160 



16 



10 



15 



20 



25 



30 



35 



and post-multiplication of a NM^. x NM^ EFFT matrix. For example, the former matrix is 
defined, 

"Y 0" 

Therefore, including the transmitter IFFT and receiver FFT operations, the input-output 
relationship is described by 



where ^ is the diagonal matrix containing the SOP bin strengths for the antenna pair (ij), 
Pre-multiplication and post-multiplication by permutation matrices Pj. and P^ represents the 
collection of all antexma combinations that correspond to a conmion frequency or SOP bin. 
This collection process, depicted in Fig. 10, results in a block diagonal matrix that relates the 
intputs and outputs: 

■H(1) 0 

P;?Y^^jGYj^^Pj. = 
which is equivalent to Equation (10). 



H(iV) 



Spatial Processing 

The spatial processing procedure is now considered. Since the SOP establishes N 
MIMO spatial chaimels that are substantially independent from one another (Equation 10), 
one can consider the spatial processing within each bin separately. Representative application 
of spatial processing to frequency bin 1 will be considered as shown in Fig. 1 1 at the 
transmitter and Fig. 12 at the receiver. Fig. 1 1 shows A/ symbols: z(l, 1) through z(l,M). 
The notation z(n,m) refers to the symbol transmitted in bin n and spatial direction m. These 
M symbols will jointly occupy frequency bin 1 . Each TSW 2 lOA-C applies a weight vector 
to the symbol appearing at its input, and the elements of the resultant vector are routed to 
Mj, summing junctions 211. One may consider the TSWs as being multipliers taking each 
input symbol and multiplying it by a vector that corresponds to a spatial direction in Mj. - 
space. Furthermore, the M vectors define a subspace in -space. Note that the TSW 
vectors are considered to be column vectors in the discussion that follows. When these M 
vectors are collected into a matrix, the result is an input orthogonalizing matrix or beneficial 
weighting matrix for that bin. For each input bm, a vector including symbols allocated to 
subchannels corresponding to the bin is multiplied by the input orthogonalizing matrix to 
obtain a result vector, elements of the result vector corresponding to the various transmitter 
antenna elements. Together, the TSWs 210A-C make up one embodiment of a Transmit 
Spatial Processor (TSP) 230. 

Each RSW 220A-C accepts Mr inputs, one from each receiver antenna path. Within 
the RSW, a weight vector is applied to the inputs (i.e. an inner-product is performed) 
thereby producing a received signal sample x(l,m): 
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x(l,m) = u(l,in) : 

where u(l,m) is the RSW for bin 1 and spatial direction m. Similar to a TSW, a RSW vector 
has an associated direction in -space. Each RSW may also be considered to be a 
multiplier. This vector is considered to be a row vector. When these M RSW vectors are 
5 collected into a matrix, the resuh is an output orthogonalizing matrix or beneficial weighting 
matrix for that bin. When a vector including symbols in a particular output bin produced by 
the SOP for each receiver anteima is multiplied by the output orthogonalizing matrix, the 
result is a vector including symbols received in that bin for various spatial directions. 
Together, the RSWs 220A-220C represent one embodiment of a Receive Spatial Processor 
10 (RSP)240. 

Through proper choice of the weight vectors applied via the TSWs and RSWs, the M 
spatial directions can be made substantially orthogonal to one another. The result is that the 
received signal sample x(l,m) depends only upon input symbols z(l,m) and not theM-l other 
input symbols for SOP bin 1. Methods for selectmg the TSP and RSP weight vectors are 

1 5 described in detail below. 

The spatial processing described above can be applied to the other N-1 fi-equency bins 
in addition to fi^uency bin 1 , The block diagram for such a system is depicted in Fig. 1 3 for 
the transmitter and Fig. 14 for the receiver. SOP processors 190 and 200 ensure that the 
fi'equency bins remain substantially orthogonal to one another while TSP 230 and RSP 240 

20 ensure that M substantially orthogonal spatial channels exist within each fi:equency bin. The 
net result is that NM substantially parallel subchannels are constructed within the MIMO 
communication system. In other words, the combination of SOP processors 190 and 200, 
TSP 230, and RSP 240 create a set of substantially independent space-frequency subchannels, 

x(/i, m) - H(n, m)z{n, m) + /(/i, m) Vn, m . 

25 This simultaneous substantial orthogonalization of space and frequency can result in a 

significant increase in spectral efficiency since multiple data streams are being communicated 
through the channel. Note that the number of substantially independent subchaimels possible, 
in the multipath case, is equal to the number of SOP bins multiplied by the number of transmit 
antennas or the number of receive antennas, whichever is smaller. Therefore, the total 

30 number of space-frequency subchannels is less than or equal toNmn (W^j. ,Mj^ ), when 
multipath is present. 

An exemplary set of TSWs and RSWs are derived from the singular value 
decomposition (SVD) of the spatial channel matrix for each bin, 

H(w) = U(/i)Z(/i)V(/i)''. 

35 The input singular matrix, V(n) , contains Mj. column vectors that define up to Mj, TSWs 
for bin w. Likewise, the output singular matrix, U(«) , contain M ^ colunrn vectors that when 
Hermitian transposed, define up to RSW row vectors for bin n. The TSWs and RSWs 
for other bins are determined in the same feshion, through an SVD decomposition of the 
spatial matrix for that bin. Using this spatial processing, substantially independent multiple 

40 streams of symbols can be transmitted and received. The strength of each subchannel is equal 
to one of the elements of the diagonal matrix S . These subchannels strengths will vary. 
Therefore, the subchannels will have varying signal to noise ratios and information capacity. 
For this reason, it may be preferable to transmit and receive only on a subset of the possible 
subchannels, or M < min{Wj. }. For example, it may be improvident to use processing 



WO 98/09385 PCT/US97/15160 



18 

complexity on the weakest subchannels that may have a very small information carrying 

capacity. In this case, spatial processing is used to increase the received power of one or 

more parallel symbol streams. It may also be preferable to use coding techniques to leverage 

strong subchannels to assist in the use of weaker subchannels. It may also be preferable to 
S allocate either bits or transmit power among the subchannels to maxinuze the amount of 

information communicated. 

The exemplary spatial processing described above requires cooperation between the 

transmitter and receiver to effectively orthogonalize the spatial channel for each bin. 

Alternatively, this orthogonalization can be accomplished at only one end of the link. This 
10 can be advantageous when one end of the link can afford more computational complexity 

than the other end. In addition, spatial orthogonalizing at one end can be advantageous when 

the channel model is known only at that end. 

Consider the case where the orthogonaIi2ation is done at the receiver. Symbols are 

transmitted along directions defined by some set of TSWs, v(n,m). When M TSWs 
15 corresponding to the same bin are collected into a matrix V(n), the composite spatial channel 

is, 

H'(n) = H(/i)V(/i). 

This composite channel describes the MIMO channel in bin n from the M-inputs to 
outputs. The spatial processing at the receiver can substantially orthogonalize this 
20 composite channel, H'(/j) , by applying appropriate RSWs even if the transmitter does no 
spatial processing. Let these RSWs be defined as the row vectors of the weighting matrix. 

Two exemplary methods for determining W^j (n) are referred to as the zero-forcing 
(ZF) solution and minimum-mean-square-error (MMSE) solution. In the ZF approach, the 
25 weighting matrix is the pseduo-(left)-inverse of the composite channel. 

This results in, 

where the identity matrix is M by M. Hence, the ZF solution, not only orthogonaUzes the 
30 spatial channel for bin n, but it equalizes the strengths of each resulting subchaimel. 

However, the signal-to-noise ratio for the various subchannels can vary widely. One skilled 
in the art will recognize that the ZF solution can result in amplification of the interference and 
noise unless the composite channel, H'(/i) , is nearly orthogonal to begin with. 

An MMSE solution, on the other hand, does not amplify noise. For the MMSE 
35 approach, the weight, W^ («) , satisfies, 

mn£}w^(«)x(«)-z(n)f}, 

or, 

(w) = R2(n)x(n)I^x(n) > 

where Rx(rt)is the covariance matrix for x(n) and R^^^^^^^jis the cross-covariance between 
40 z(n) and x(n). Using, 

x(n) = H'(/i)z(/2) + r(/i), 
and the fact that R,^„j = a^l .results in the MMSE weight. 



f 1 
W^ (/?) = H'(n)'' lI\n)RXnf + -yR,. 
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Note that when I(n) is spatially white noise, then Rj = cT; I . 

Similarly to the above orthogonalization at the receiver, the channel can be orthogonalized at 
the transmit end only. For this to occur, the transmitter is required to have knowledge of the 
5 RSWs to be used by the receiver. In a TDD channel, where the channel exhibits reciprocity, 
these RSWs can be learned -when that transceiver uses TSW directions equal to the RSW 
directions. Alternatively, the receiver may not do any spatial processing, so the transmitter is 
responsible for spatial orthogonalization. 

10 In this case, the composite channel is 

H'(n) = U(«)H(n), 

where the matrix U(n) is composed of the RSW row vectors, u(n,m). This composite 
channel describes the MEMO channel in bin w from A/j. inputs to M outputs. Similar to the 
previous case, the transmitter can substantially orthogonalize this composite channel, H'(«) , 
1 5 by applying appropriate TSWs. These TSWs are the colunm vectors of the weighting matrix, 

The transmit weighting can be determined using the ZF or the MMSE approach. In 
the ZF approach, the weighting matrix is equal to the pseudo-(right)-inverse of H'(«) . The 
MMSE solution satisfies 
20 mm £iH'(H)W,(/i)z(/i) + lX«)-z(n^ |f ), 



An important simplification to the general space-fi-equency processing technique is the use of 
only one spatial direction for each bin of the SOP. This case is depicted in Fig. 1 S for the 

25 transmitter and Fig. 1 6 for the receiver. In this case, only N subchannels are created. The N 
input symbols, z(l,l) through z(N,l), are processed by N TSWs 210A-B that weight and 
allocate these N symbols among the A/j. identical SOP processors 190. At the receiver, the 
antenna samples are processed by SOP processors 200. TheM^j SOP outputs 
corresponding to a common bin are weighted and combined in N RSWs 220A-B. With N 

30 such weightings, the result is N outputs, x(l, 1) through x(N, 1), of the N substantially 
orthogonal subchannels. 

When only one spatial direction is used in the TSP and RSP, one exemplary choice for 
the particular weightings are the TSW and RSW directions that result in maximum 
subchannel strength. This maximizes the signai-to-noise ratio (S>]R) of the received signals, 

35 x(l,l) through x(N, 1). In this case, the optimal weightings should satisfy, 

max |u(/i)H(/i)v(w)|', 

U(||).T(||)' ' 

with the implicit constraint on the size (2-norm) of both the RSW weight u(n) and the TSW 
weight v(n). To determine the solution to this optimization problem, consider the SOP 
outputs for bin n when a single TSW, v(n,l), is used, 

40 : = H(rt)v(w,l)z(/f,l) = h(»)z(/2,l) . 
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The quantity h(n) is referred to as the received channel vector. A channel 
identification technique is used to determine the received channel vector. Therefore, the 
optimal RSW weight is equal to the Hermitian of the received channel vector, h(n), 

u(/j) = h(/i)'^. ; 
5 Note that this is true regardless of the particular value of v(n). The optimal TSW direction, 
on the other hand, satisfies, 

max u(/j)H(«)v(n) = max v(n)^ K(n)^ H(w)v(n) , 

where the optimal RSW direction has been used. The optimal TSW for bin n is equal to the 
scaled maximum eigenvector of the matrix B,{n)^ H(/7) . One skilled in the art will recognize 

1 0 that the optimal RSW is also equal to the scaled maximum eigenvector of H(«)H(n)^ . 

A fiirther advantageous simplification of the above techniques is the use of one or 
more common TSW and RSW directions for all bins. In other words, every bin has the same 
TSW and RSW weights. These weight vectors may also consider delay elements. In one 
embodiment, these weights are determined to maxunize the SNR of the received signals, 

15 averaged over fi^equency «. This is depicted in Fig. 17 for the transmitter and Fig. 18 for the 
receiver. Consider this embodiment with one spatial direction. In this case, the TSW and 
RSW weights satisfy, 

max£jtfH(»)vf , 50 

tl,V 

Note that the expectation operator, E„ , represents averaging over SOP bins. TMs averaging 
could also be done over multiple bursts in addition to fi-equency. The solution to this problem 
is when v is equal to the maximum eigenvector of 

R,^E„^^in)H{n)l 51 

and u is equal to the maximum eigenvector of the covariance matrix formed firom averaging 
the outer product of the receive vector channel, 
25 R,=£:>(«)h(/i)^}. 52 

The quantity R^is the spatial covariance matrix that describes preferable directions to 

transmit to the desired receiver, a desired subspace. 

This technique can be generalized to the case where multiple directions are utilized. 

In this case, MTSWs and MRS Ws are determined to maxunize the average (over bin) SNR 
30 received through the M spatial directions. The M spatial du-ections will not necessarily be 

orthogonal to each other. Therefore, there will be spatial crosstalk m the received symbols. 

Multidimensional encoding and decoding techniques discussed below can then achieve a 

multiplicative rate increase in the presence of such crosstalk. 

Alternatively, the receiver can spatially orthogonalize the subchannels by fiirther 
35 weighting of the M outputs fi"om the RSWs. The composite spatial channel at bin n, with the 

RSWs and TSWs included is 

H'(/i) = UH(w)V, 53 
where the matrix U is made up rows equal to the RSW directions and V is a matrix with 
columns equal to the TSW directions. Since U and V were determined based on an average 
40 SNR criterion for all bins, the composite matrix H'(«) will not be diagonal. Hence, the 

receiver can apply the additional weight, W^ (n) , to orthogonalize H'(«) . Alternatively, the 
transmitter can use the additional weight, Wj. («) , to spatial orthogonalize the composite 
channel. Exemplary solutions for these weightings are the joint SVD, the ZF and MMSE. 
The advantage of this approach is that the processing required to adapt all N SOP bin matrix 
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channels may be substantially higher than the processing complexity to adapt the average TSP 
andRSP. 

The rejection and prevention of interference can be accomplished in conjimction with 
the space-frequency processing discussed above. This is especially useful when the numbo* of 
S spatial directions used for communication is less than the number of antennas. This case 
occurs when weak spatial directions are not utilized or when the number of antennas at the 
receiver and transmitter are not the same. In either case, one or both ends of the 
communication link have extra spatial degrees of freedom to use for the purpose of mitigating 
mterference. 

1 0 The amount of interference arriving at an antenna array can be quantified by the 

interference covariance matrix, 

where I{n) is the length interference plus noise vector received in SOP bin n. This 
matrix defines an imdesired interference plus noise subspace in -space for bin n. The 

1 5 interference plus noise energy that contaminates a particular received subchannel symbol with 
bin n and spatial direction m, is equal to, 

u^(w,w)Ri(/i)u(/i,ot), 
where u(n, m) is the combining weight vector for the RSP(n,m). An advantageous 
interference rejection technique is then to ^'whiten" the effect of the interference across the 

20 spatial directions, so that the interference is minimized and spread evenly across all spatial 
directions used. Therefore, each of the RSP weighting vectors are modified by the matrix 

^ r;''^(/i), 

u'(/i,m) = u(w,w)R-'''(/i). 
Alternatively, the RSP weighting vectors are the vectors of the output singular 
25 matrix, JJ'Qi) , from the SVD of the modified spatial channel, 

R-'^" (w)H(«) = J]Xn)r(n)VXnf . 
Note that a very usefiil simplification of the above interference rejection technique is to 
average the interference covariance matrix over all N bins and possibly a set of bursts to 
arrive at an average spatial interference covariance matrix, R^ , that is independent of bin n. 
30 In this case, every RSP combining vector is modified in the same way due to interference. 
This approach can significantly reduce the amount of computations needed to determine R^ 

and RJ^^^ . Note that it is often beneficial to add a scaled matrix identity term to estimates of 
the interference covariance matrix to reduce the sensitivity of these interference mitigation 
approaches to covariance estimation errors. 

35 Similar interference mitigation techniques can be advantageously employed at the 

transmitter to reduce the amount of interference radiated to xmintentioxud receivers. In the 
TDD channel, reciprocity in the radio link allows the undesired receive interference subspace 
in each SOP bin to be accurately used to describe the transmitter subspace. That is, the 
amount of interference transmitted to unintentional receivers is 

40 (/I, m)Rj («) v(n, m) , 60 

where v(n, m) is the transmit weight vector for the TSW(n,m). An optimal interference 
reduction approach is then to minimize and *Svhiten" the transmitted interference across 
spatial directions. In the same fashion as the receiver case, the TSW vectors are modified by 
the matrix Ry*'^ (w) . Alternatively, the TSP weight vectors are the vectors of the input 

45 singular matrix V'(rt) from tiie SVD of the modified spatial channel. 
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n{n)R-''\n) = V\n)lXfi)yXnf . 61 
Again, a significant simplification occurs when the interference covariance matrix is 
determined by averaging over frequency or SOP bins. It is especially advantageous to 
average over SOP bins in a frequency-division-duplex (FDD) system, where significant 
5 averaging of the receive convariance matrix results in a good estimate of the transmit 
covariance, even though instantaneous channel reciprocity does not hold. 

Interference rejection at the receiver and interference reduction at the transmitter are 
done together by simply combining the two techniques outlined above. In this case, the RSP 
vectors and TSP vectors are contained in the input and output matrices of the SVD o^ 

10 R7i;'(/»)H(/i)R7y^(«), 

As outlined previously, it can be advantageous to use the same TSWs and RSWs for 
. all bins. This approach can be combined with interference mitigation by the determining the 
transmit and receive weight vectors that maximize average power delivered to the receiver of 
1 5 interest, while at the same time, minimizing power delivered to other imdesired receivers. 
There are various optimization problems that can be posed to determine these TSP or RSP 
directions, each involving the desired receiver covariance matrix and the undesired covariance 
matrix. For example, one TSP problem is 

max v^R^v such that v"R^ v ^ Pj and v^v <. Pj. . 

20 That is, a TSP direction is chosen for all SOP bins that transmits the maximum amount of 
r-..,. power to the desired receiver while maintaining a transmit power limit, Pj> , and a transmitted 
interference limit, Pj . For this particular problem, the TSP direction is equal to the maximum 
generalized eigenvector of the matrix pair { , (Rj /P^ + 1 /P,. ) }. One example of an 
eflfecdve interference rejecting RSP for all SOP bins is a weighing that maximizes the average 
25 received SINR. The RSP that maximizes SINR is the maximum eigenvector of the matrix 

r;'"r,R7"V. 

One fiirther simplification to the above algorithm is to model the interference and/or 
the desired covariance as diagonal, or neariy diagonal. When both R^ and R^ are diagonal, 
the solution to the above optimization problem reduces to the maximal ratio SINR combiner, 

30 u, and transmitter, v. It is also sometimes preferable to only consider other subsets of the 
elements of either the desired or interference covariance matrices. 

One skilled in the art will also recognize that all the TSPs and RSPs can be used when 
there is only one SOP bin, such as a common frequency. 
Space Frequency Coding 

35 Many of tiie advantageous space frequency encoding techniques embodied in the 

invention may be broadly classified in two exemplary categories. The first category involves 
techniques wherein the spatial matrix channel within each SOP bin undergoes space frequency 
processing at the transmitter, or the receiver or both, resulting in a substantially independent 
set of one or more parallel communication subchannels within each SOP bin. The objective 

40 of the encoder and decoder in this case is to appropriately allocate the transmitted information 
among multiple independent space-frequency subchannels using interleaving, power and bit 
loading, or both. The second category of space frequency encoding involves transmitting and 
receiving one or more symbol sequences in each SOP bin using one or more transmitter and 
receiver weight vector combinations that are not necessarily intended to create independent 

45 spatial subchannels within each SOP bin. This results in significant cross-talk between each 
symbol stream present at the receiver output. A decoder then uses knowledge of the 
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equivalent matrix channel within each SOP bin, and knowledge of the set of possible encoder 
sequences to estimate the encoder symbol sequence that gave rise to the cross-talk rich 
output SOP bin vector sequence. The main diflFerentiating feature between the first approach 
and the second approach is the presence or lack of spatial processing that results in 
5 substantially orthogonal spatial subchannels within each SOP bin. Both approaches have the 
advantageous ability to multiply the data rate that can be achieved in MIMO channels with 
multipath. 

Coding for Substantially Orthogonal Space-Frequency Subchannels 

10 In applications where the spatial channels are processed to achieve multiple 

substantially independent spatial subchannels within each SOP bin, an advantageous 
embodiment of the invention involves encoding the input data sequence into a digital symbol 
stream that is then routed in various beneficial ways through the available parallel space 
fi-equency subchannels. Fig. 22 depicts a preferred embodiment. This embodiment involves 

1 5 distributing the symbol outputs of a single encoder among all of the available space firequency 
subchaimels. Several known coding schemes that can be combined effectively with space 
fi'equency processing to distribute information transmission over the space and fi-equency 
dimensions of a communication channel This discussion assumes estimation of the MIMO 
channel by transmitting a series of training symbol sequences fi*om each antexma element as 

20 discussed herein. The discussion fiuther assumes that the receiver and transmitter either 
share the information required to decompose the channel into parallel sub-channels, or the 
TDD techniques discussed herein are used to do the same. 
■f Referring again to Fig. 22, the preferred embodiment exploits a three layer coding 

system. The first layer of codmg includes the combination of transmitter TSWs 2 IDA 

25 through 210B, Transmitter SOP processors 190, receiver SOP processors 200, and receiver 
RSWs 220A through 220B. This first layer of coding performs the spatial processmg. The 
second layer of codmg includes a Trellis Encoder and Inteileaver (420) at the transmitter in 
combination with a Deinterleaver and ML Detector 430 at the receiver. The third layer code 
involves Reed Solomon (RS) Encoder 410 at the transmitter in combination with an RS 

30 Decoder 440 at the receiver. The bit level RS coding occurs prior to the trellis encoding and 
the Reed Solomon codeword detector acts upon the bit sequence fi'om the ML detector. The 
fourth layer of coding involves an ARQ code that recognizes Reed Solomon codeword errors 
at the receiver in the Receiver ARQ Buffer Control 450 and requests a codeword 
retransmission fi'om the Transmitter ARQ Buffer Control 400. The retransmission request is 

35 made through a Reverse Link Control Channel 460. The reverse control channel is a well 
known radio system concept and will not be discussed herein. This combination of coding 
techniques and space fi-equency processing is preferable because it provides for a rich 
combination of space and fi-equency diversity and it is capable of obtaining very low bit error 
rates. The detailed operation of the RS encoder and decoder, as well as the ARQ system is 

40 well known to one sldlled in the art. Following this discussion, it will be clear to one skilled 
in the art that other combinations of one or more of these four coding elements may be 
employed with advantageous results in various applications. 

The trellis coding step may be substituted with CBM-QAM or a turbo code. 
Similarly, the Reed-Solomon code may be substituted with a block code, or with an error 

45 checking code such as a CRC code. The transmitter end would then include the necessary 
encoder and the receiver end would include the necessary decoder. 

There are at least two basic methods for employing trellis coding to distribute 
information among substantially independent space fi'equency subchannels. One method is 
adaptive encoding that modifies the bit and power loading for each subchanel according to its 
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quality. The second method involves maintaining constant power and bit loading for all space 
frequency subchannels. Both ofthese methods are discussed below. 

Space Frequency Trellis Coding with Orthogonal Spatial Subchannels and Adaptive 
5 Power and Bit Loading 

Fig. 22 depicts the coding and interleaving detail for the transmitter and receiver 
portions of the present embodiment. Encoding and Interleaving system 10 encodes the data 
into a set of complex symbols. Each of the complex symbols is then allocated to a particular 

10 transmitter TSW (210A through 210B). The input to each TSW forms a vector of frequency 
domain symbols that are fed into the same bin of one or more transmit SOP processors. Each 
transmitter TSW, possibly in conjunction with a receiver RSW converts the matrix channel 
within each SOP bin into a set of substantially orthogonal space frequency subchannels using 
one of the methods discussed herein. 

15 Fig. 23 displays a more detailed diagram of the encoder and interleaver. An 

Information Allocation Unit 360 assigns the bits and. the transmitter power that will be 
allocated to each space frequency subchannel. One method for accomplishing this assignment 
is the so-called gap analysis. In this technique, a particular coset code with an associated 
lattice structure is characterized by first determining the SNR required to achieve a theoretical 

20 capacity equal to the desked data rate. The code gap is then the SNR multiplier required to 
achieve the target probability of error at the desired data rate. In a parallel channel 
communication system, this gap can be used to determine the power and bit distributions that 

^ maximize data rate subject to a probability of error constraint. With a coding gap of a , the 
rate maxunizing water-filling solution for the space fi^quency subchannels becomes 



25 piji^m)- 



f 2 ^ 



\Xin,mf) ' 

where is the noise power and m is the spatial index and n is the DFT frequency index. 
The bit allocation per sub-channel is then given by 

After the power and bit loading assignments are accomplished in the Information Allocation 
30 Unit, the bits are encoded with a Trellis Encoder 370. It is not possible to achieve infinhe bit 
resolution (granularity) with coset codes. Therefore the gap analysis solution should be 
modified. Several bit loadmg algorithms exist to resolve this problem. One method mvolves 
rounding down the water filling solution to the nearest available quantization. The granularity 
of possible bit allocations is determined by the dimensionality of the coset code lattice 
35 structure. In the MIMO channel communication structures described herein, the orthogonal 
constellation dimensions are the complex plane, space, and frequency. 

Fig. 26 illustrates an example of a practical method for bit loading with a trellis 
encoder that uses a one dimensional QAM symbol constellation. The bit load is adjusted 
down for a ^ven trellis encoder output symbol by assigning a number of fixed zeros to one or 
40 more of the input bhs to the encoder. Fig. 26 shows the operation of the trellis encoder and 
the trellis state diagram for the decoder for four successive symbol transmissions. Four bits 
are assigned to the first space-frequency subchannel. A first subchannel trellis encoder input 
350 is assigned 4 bits so Symbol 1 can take on any one of 32 values. There are two bits 
feeding the convolutional encoder, and two bits feeding the coset select. The ML detector at 
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the receiver uses the trellis state diagram and the channel state information to solve the 
maximum likelihood recursion. This is eflBciently accomplished with the \^terbi algorithm. 
The trellis code state diagram defines a set of symbol sequence possibilities {Z}. The space 
fi-equency subchannel is denoted H{n, k) , for SOP bin n at burst k. The maximimi likelihood 
5 equation is then given by 

where z(n) is the symbol hypothesized for SOP bin w. The decoder output state diagram for 
first space fi-equency subchannel 340 includes four possible parallel transitions for each trellis 
branch and all of the trellis branches are possible. The second space fi-equency subchannel in 

10 the sequence is assigned three bits so a second trellis encoder input 352 shows one bit fed 
into the coset select with two bits still feeding the convolutional encoder. A decoder state 
diagram 342 for the second space firequency sufachaanel has only two parallel transitions for 
each trellis branch but still maintains ail trellis branch possibilities. Continuing in succession, 
a third space fi-equency subchannel is assigned only two bits to an encoder input 354 so there 

15 are no parallel transitions considered by a trellis decoder state diagram 344. In a fourth space 
fi-equency subchannel, only one bit is assigned to an encoder input 356 so there are no parallel 
transitions and some of the trellis state branches (346) are no longer considered by the 
decoder. It is understood that Fig. 26 is provided as a graphical aid and is not intended to 
represent an actual design. 
^20 The maximimz Euclidean distance error sequence design metric is one preferable 

choice for a trellis encoder used with the parallel space firequency chatmel with this bit and 
power loading embodiment of the invention. Other code error sequence design metrics that 
are advantageous in various application conditions include product distance and periodic 
product distance. 

25 Referring again to Fig. 23, the output of the encoder is interleaved across the various 

space fi-equency subchannels using Interleaving block 260. Typically, the interleaving process 
cUstributes the symbols so that symbols that are near one another at the encoder output are 
well separated in both the SOP bin assignment and the spatial subchannel assignment. This 
distributes the effects of channel estimation errors and localized fi'equency domain or spatial 

30 domain interference so that the decoder error is reduced. It is understood that the bit and 
power assignments by Information Allocation block 360 take place with knowledge of the 
post-interleaved chamiel strength. It is understood that the encoding and decoding process 
can begin and end within one burst, or it may take place over a multitude of bursts. 

One skilled m the art will recognize tfiat a multitude of less sophisticated adaptive 

35 power and bit loading algorithms can be advantageously applied to a substantially 

independent set of space fi'equency subchannels. One example is an algorithm wherein a 
space fi-equency subchannel is either loaded with maximum power or no power and the bit 
distribution may be adjusted in only two increments. 

A second alternative embodiment shown m Fig. 19 includes one encoder for each 

40 SOP bin, with the output symbols of each encoder allocated among several spatial 

subchannels. A third embodiment shown in Fig. 20 involves one encoder for each spatial 
subchannel, with the output symbols of each encoder distributed among the SOP bins for that 
spatial subchannel. A fourth embodiment shown in Fig. 21 involves a separate encoder for 
each available space fi'equency subchannel. 

45 It will be clear to one skilled m the art that the channel estimation tools taught herein 

are very usefiil in improving the accuracy of the channel estimates used for the bit loading and 
decoding process. 
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One skilled in the art will recognize that many of the other coding techniques for 
parallel sub-channel bit loading communication systems, not mentioned here, can also be 
applied to the present invention. 

S Space Frequency Trellis Coding with Orthogonal Spatial Subchannels and Flat Power 
and Bit Distribution 

In some cases it is difiScuh to adaptively load the power and bit assignments for each 
available space frequency subchannel. For example, the transmitter and receiver may not be 

1 0 able to adapt the loading fast enough to accommodate time domain variation in the channel. 
In another example, the required feedback from the receiver to the transmitter requires a 
significant portion of the available reverse link bit rate. Adaptive bit loading may also be 
overly compUcated for certain applications. Thus, it is often advantageous to encode and 
decode a symbol stream in such a manner that the power and bit allocation is constant for all 

1 5 space frequency subchannels. This is easily accomplished by employing the embodiments 
depicted in Figs. 22-23 , and assigmng a constant power and bit allocation to all space 
frequency subchannels in the Information Allocation block 360. 



20 Space Frequency Coding Without Orthogonal Spatial Subchannels 

In applications where the spatial channels are not processed to achieve substantially 
orthogonal spatial subchannels v^thin each SOP bin, an advantageous embodiment of the 
*>- invention involves utilization of a vector maximum likelihood decoder in the receiver to 

decode a symbol sequence that includes multiple symbols per SOP bin. The vector maximiun 
25 likelihood detector is capable of determining tfie transmitted symbol vector in each SOP bin 
even in the presence of spatial subchannels that contain significant cross-coupling between the 
chaimels. The vector maximimi likelihood detector uses an estimate of the matrix channel 
from each SOP bin to decode a sequence of groups of symbols with one group for each SOP 
bin. The groupings will be referred to here as a multidunensional symbol vector, or simply a 
30 symbol vector. The ML detector uses an estimate of the matrix channel that exists in each 
SOP bin to find the most likely sequence of transmitted encoder vector symbols. 

Fig. 24 depicts a transmitter system wherein multiple space/frequency subchannels are 
employed without spatial orthogonalization. Fig. 25 depicts a receiver sjrstem for this 
application. 

35 The bit sequence b(k) is encoded into a sequence of multidimensional symbol vectors 

in a Bit to Symbol Encoding block 250. Each output of the encoder is an Mo by 1 complex 
symbol vector, where Mo is the number of spatial directions that will be used for transmission. 
Note that Mo is preferably chosen to be less than or equal to Mt. A preferable construction 
of the encoder is a multidimensional trellis encoder. One advantageous metric for designing 

40 the trellis encoder constellation and convolutional encoder polynomial will be provided 
below. Within the previously discussed Symbol Interieaver block 260, the vector symbol 
sequence is demultiplexed and interleaved with a Symbol Sequence Demultiplexor 300 and a 
Transmit Symbol Routing block 310. Transmit Symbol Routing block 310 interleaves the 
vector symbol sequence so that the elements of a given vector symbol are grouped together 

45 and transmitted in one SOP bin. Thus, diflFerent vectors are separated by a multitude of SOP 
bins before transmission, but all elements within the vector symbol share the same SOP bin. 
The purpose of the interieaver is to distribute the vector symbol sequence so that the fading 
present in the matrix channels within the SOP bins is randomized at the output of the receiver 
interieaver. The decoder can recover information associated with symbols that are 
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transmitted through SOP bins that experience a deep &de, provided that the adjacent symbols 
do not also experience the same fade. Since there is often a high degree of correlation in the 
fading experienced by adjacent SOP bins, the interieaver makes the fading more random and 
improves decoder error performance. After interleaving, each element of a vector symbol is 
5 assigned to one antenna for the SOP bin assigned to that vector symbol. Transmitter SOP 
processors 190 perform the transmitter portion of the SOP. 

After the transmitter SOP, it is often advantageous to perform spatial processing with 
TSP 230. It is understood that the matrix representing the operation of TSP 230, i.e., the 
Transmitter Weight Matrix, may also be an identity matrix so that no weighting is 

10 implemented. It can be beneficial to choose a number of spatial directions that is less than the 
number of transmitter antennas. In this case, the Transmitter Weight Matrix increases the 
dimensionality of the time domain vector sequence from the SOP bank. As an example of 
when it is advantageous to choose a subset of the available transmitter spatial directions, if 
the receiver has fewer antennas than the transmitter, then it is known that the information 

1 5 capacity of the matrix channels within each SOP bin will not support a number of parallel 

information subchannels that is greater than the number of receive antennas. This implies that 
the number of symbols in each transmitted symbol vector, and hence the number of 
transmitted spatial directions, should not be greater than the number of receiver antennas. As 
another example, in a Rayleigh fading channel, the smallest singular values of an Mr by Mr 

20 matrix channel are on average much weaker than the largest singular value. This implies that 
the average information capacity contained in the smallest singular value may not justify the 
extra signal processing complexity required to transmit over that dimension. In both of these 
,4 cases, it is advisable to choose an advantageous subset of the available transmit spatial 
directions. 

25 The transmitter may not have knowledge of the individual channel matrices within 

each SOP bin but may have knowledge of the covariance statistics of the channel matrices, 
averaged over frequency, or time, or both. In such cases, the Transmitter Weight Matrix can 
be optimized to select one or more spatial directions that maximize the average received 
power for the chosen number of spatial directions. The procedure for optimizing the 

30 Transmitter Weight Matrix for this criteria is defined by Equations 50 to 52 and the 

associated discussion. This is one preferred method of selecting an advantageous set of 
spatial directions for the Transmitter Weight Matrix. Another advantageous criteria for 
selecting the transmitter spatial directions is to maximize average received power subject to 
constraints on the average interference power radiated to unintentional receivers. The 

3 5 procedure for optimizmg the Transmitter Weight matrix for this criterion is defined by 
Equations 60 and 61 and the associated discussion. 

After the time domain signal is spatially processed, the signal is upconverted to the RF 
carrier frequency using Modulation and RF System 40 before being radiated by Transmit 
Antennas 5 1 . Referring now to Fig. 25, at the receiver the signal is downconverted and 

40 (hgitized by Receive Antennas 1 1 1 and Demodulation and RF System blocks 120. The RSP 
240 may then be used to process the time domain signal. The operation of RSP 240 may be 
characterized by a Receiver Weight Matrix which may be an identity matrix. One 
embodiment involves optimizing the RSP weights to reduce the number of received signals 
from Mr to Mo, which is the ninnber of elements in the transmitted symbol vector and is also 

45 the number of transmitted spatial directions. In this case, the Receiver Weight Matrix can be 
optimized to increase the average signal power in each received spatial direction. The 
optimization procedure to accomplish this is defined by Equations 50 to 53 and the associated 
text. 
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Channel ID block 130 is used to estimate the matrix channel in each SOP bin. 
Procedures for channel estimation are described below. Channel state information for each 
SOP bin is fed into a Symbol to Bit Detector 280 which decodes the symbol sequence after it 
is passed through a Symbol Deinterleaver 270. 
S At the receiver, after de-interleaving the SOP bins, the space-fi'equency sequence is 

again converted into a serial symbol stream by Demultiplexor 300. For a given set of spatio- 
temporal vector symbol sequraice posabQities {Z}, and an estimate, H(n, A) , of the channel 
matrix in each SOP bin n at burst k, the maximum likelihood detector is given by equation 
(70): 

10 ^(l)^z(2)^..^z(A^)^}=arg^^^^^^^^^^^^ 

where z(n) is the vector representing the code segment hypothesized for SOP bin «, and 

Rj(n,k) is the estimated noise plus interference covariance matrix for SOP bin n and time k. 

This equation can be solved eflBciently using a vector ML detector. The SOP fain channel 

matrix estimates are understood to include the effects of the Transmitter Weight Matrix and 
1 5 the Receiver Weight Matrix. It is understood that the noise pre-whitening step in the ML 

detector cost ftmction can be substituted by a bank of RSPs that perform the interference 

whitening as described herein. 

In a Rayleigh fading channel, a desirable metric for designing the trellis code is given 

by the product of a sum involving the two-norm of vector segments of the trellis code error 
20 sequence: 

where q is the number of SOP bins in the error sequence, and e(n) is the vector difference 
between the true multi-dimensional code symbol segment and the incorrect multi-dimensional 
symbol code segment for SOP bin w. This code design metric is a generalization on the 

25 conventional product distance metric which contains a scalar error entry in the product 

equation while the new code design metric contains a vector two norm entry in the product 
equation. It should now be evident that the multidimensional encoder can be realized by 
either directly producing a vector consisting of a multidimensional QAM symbol with the 
encoder output or by grouping complex QAM symbols from a one dimensional encoder 

30 output into a vector. The vector symbol encoder alternative is preferred in some cases 
because this approach provides for a larger metric search result and hence a better fading 
code. After deinterleaving, the decoder that is used at the receiver searches over all possible 
multidimensional symbols within each SOP bin to maximize Equation 70. It is understood 
that one skilled in the art will recognize after this discussion that other desirable metrics such 

35 as Euclidean distance metrics, metrics designed for Ridan fading channels, periodic product 
distance metrics, and others are straightforward to construct and space-frequency codes can 
then be determined through well known exhaustive search techniques. 

In either the one dimensional encoder case, or the multidimensional encoder case, the 
encoder constellation selection and code polynomial search to maximize the metric can be 

40 carried out using a number of well known procedures. 

It is possible to improve the performance of the space-frequency coding system 
described above by using a number of transmitter antennas, or a number of receiver antennas, 
that is greater than the number of symbols transmitted in each SOP bin. If the number of 
receiver antennas is greater than the number of symbols in each SOP bin, then simply applying 

45 the approach described above is advantageous. If the number of transmitter antennas is 
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greater than the number of symbols transmitted in each SOP bin, then the techniques 
embodied in Equation 70 are advantageous. 

Channel Identification 

5 The operation of Channel Identification block and Training Symbol Injection block 20 

will now be described. The transceiver should determine the MIMO channel in order to form 
the TSWs and RSWs. For coherent spatial processing and detection, the receiver should 
obtain an estimate of the channel. We wish to identify the set of matrix channels that results 
after processing by the transmitter and receiver portions of an SOP. The notation for this 

10 channel is H(/7) V n where n is the SOP bin index. Channel identification techniques 

embodied herein can be applied to several preferable SOP pairs including the IFFT-FFT with 
cyclic prefix, the multiband filter bank, or any other of a number of well-known SOPs. The 
following exemplary channel identification approach exploits the correlated fi:'equency &ding 
across and possibly the correlated time fading in the channel. The correlation in the fi:^uency 

15 domain arises due to the limited time delay spread of the multipath channel. The correlation 
in time is due to the fact that the channel, while tune-varying, is driven by band-limited 
Doppler firequencies created by objects, which can include the transmitter and/or receiver, 
moving in the physical environment. 

The wireless link is bidu-ectional, therefore each end of the link should estimate not 

20 only a receive channel, but also a transmit channel. For example, a base station should 
estimate both an uplink and downlink channel. In systems which employ time division 
-^K duplexing (TDD), electromagnetic reciprocity implies the receive and transmit propagation 
environments are the same, allowing the transmit channel to be estimated firom the receive 
channel. However, the transmit and receive electronic responses are not necessarily 

25 reciprocal, and because the net channel response includes the electronics, a calibration 
procedure should be used to account for these differences. This calibration procediae 
provides for matching in the amplitude and phase response between the multiple transmitter 
and receiver fi-equency converters. Several TDD calibration procedures are known in the 
prior art and will not be discussed herein. 

30 In systems employing firequency division duplexing (FDD), the propagation medium is 

not reciprocal; however, the paths* angles and average strengths are the same for transmit and 
receive. This enables the use of subspace reciprocity, but incurs a more rigorous calibration 
requirement. The FDD calibration should insure subspace reciprocity which requires that the 
array response vector at a given angle on receive is proportional to the corresponding vector 

35 on transmit. This requirement is satisfied by again calibrating the amplitude and phase 
differences among the multiple transmit and receive frequency converter channds and by 
matching the transmit and receive antenna element response as well as the array geometry. 

An alternative approach to transmit channel estimation in FDD systems uses feedback. 
The transmit channel is measured by sendmg training symbols to the receiver, which records 

40 the amplitude induced by the training symbols. Using receiver to transmitter feedback on a 
separate feedback control channel, the training responses are sent back to the transmitter. The 
transmitter, knowing the trainmg excitations it used and the corresponding responses through 
feeedback, the forward channel can be estimated. 

In general, channel identification can be done either with or without training. A 

45 desirable channel identification algorithm should be robust to and operate in a variety of 

modem implementations. A preferable MIMO channel identification technique operates with 
embedded training inserted into the data stream by Training Symbol Injection block 20 in 
each burst. In this case, both data symbols and training symbols may be transmitted within a 
single burst. Furthermore, the channel can be determined in one burst, or filtering training 
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data gathered over multiple bursts. Being able to update the channel estimates after every 
received burst makes the overall communication system robust to time variation in the 
channel In addition, frequent channel estimates reduce the destructive effects of imperfect 
carrier frequency recovery. Since imperfect carrier recovery imparts a phase shift to the 
S charmel that continues to grow with time, shortening the time between charmel estimation 
events keeps the channel estimation information from becoming "stale". Note, however, any 
of the well known blind chaimel estimation techniques can be used to determine the training 
symbol outputs as an alternative to using training. However, adaptive blind training is more 
prone to generating burst errors, 
10 The parameters to be identified are the N MIMO spatial charmel matrices. Hence, 

there are N'M^ -M^ complex elements to be determined, 

V;;€[l,;\^],V/G[l,Mj,V7e[l,Mr] 

By exploiting whatever correlation exists across the SOP bins, it may be possible to reduce 

the amount of overhead required to identify the channel. The amount of correlation that 
1 5 exists across SOP bins is determined by the specific implementation of the SOP. If the SOP 

implementation includes the IFFT-FFT pair, and the length of the FIR channel is time limited 

with V « AT, then a relatively large degree of correlation exists across the SOP bins. 

In certain embodiments of the invention, the desired technique should identify the 

MIMO channel on a burst-by-burst basis, such as those with rapidly time-varying channels. 
20 This implies that training data should be included in every burst. If the throughput of 

information is to be maximized, the amount of training data in each burst should be 
V ^: ^5 minimized. It is therefore useful to determine the minimum amount of training data required, 

per burst, that allows full characterization of the channel by the receiver. It turns out that the 

minimum number of training symbols required to suflSciently excite the MIMO chanriel for 
25 estimation with an OFDM SOP is M^v . To understand this result, consider the identification 

of a SISO charmel, where each of the N values of the vector H, y should be found. These N 

values are not independent since, 

Hw=y[J] = X0-'Z, 

where X is a vector of all SOP bm outputs for anterma / , Z is a vector of bin inputs for 
30 antenna 7, and h is a vector of the time-domain FIR channel from anteima j to antenna / . 

The matrix operator ®"* represents element by element divide. Since the time-domain 
channel is time limited to v samples, only v values of the transmitted symbols, Z , need to 
be training values. Furthermore, the identification of the SIMO channel only requires the 
same set of v transmitted traming tones, since each SISO component in the SIMO channel is 
35 excited by the same input data. In a system embodiment with multiple inputs (A/j. > l), 
identification of the MIMO channel requires the identification of Mj separate SIMO 
channels. Hence, only M^v training symbols are needed to suflSciently excite the MIMO 
charmel for chaimel identification. 



40 MIMO identification 

The identification of the MIMO chaimel is accomplished by separately exciting each 
of the transmit antermas that will be used for communication. This decomposes the MIMO 
identification problem into Mj. SIMO identification problems. In order to accomplish 
channel identification in a single burst, M^. mutually exclusive sets of vbins are selected from 
45 the N available bins to carry training symbols. Each transmitter anterma carries training 
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symbols in a unique one of the Mj, sets of bins, while transmitting no energy in the bins 
contained in the imion of the remaining A/j. - 1 sets of v bins. This is accomplished by 
choosing the TSWs 210A-C that correspond to training bins such that a single entry in the 
vector is "1" and the remaining entries equal to "0". It is the entry of a TSW that is set 
5 equal to "1" for those training symbols which are to be transmitted from the antenna. For 
example, say that symbol bin ;? = 2 is one of the training bins ^sociated with transmit 
antennas. Then, 

TSW(2,1)=[0 0 1 0 - 0r,andTSW(2,m) = 0for Vm;tl, 
and the corresponding training symbol 2(2, 1 ). By examining the contents of each set of 
1 0 training bins separately, the MIMO channel response is determined by finding Mj. 
independent SIMO channel responses. 

In embodiments in which rapid updates of the channel estimate is not required, 
another exemplary training scheme may be employied. This training scheme involves using 
just one set of v training bins. On a given burst, one of the transmit antennas sends training 
1 5 symbols in the training bins and the other antennas transmit no energy in those bins. This 
allows the receiver to identify one of the Mj SIMO channels. On the next burst, a dififerent 
antenna sends training s}rmbois in the trainmg bins while the other antennas transmit no 
energy in those same bins. The receiver is then able to idmtify another set of N SIMO 
channels. This procedure is repeated until training data has been sent by each of the transmit 
20 antennas, allowing the entire MIMO channel to be identified. The entire procedure is 
repeated continuously so that fiill channel is determined every Mj. bursts. 



SIMO channel identification 



25 We have just shown that identification of the MIMO channel can be accomplished by 

successive identification of each SIMO channel. It is therefore usefiil to discuss specific 
techniques for obtaining a SIMO channel response. The foUowiiig discussion assumes that 
the SOP is the IFFT-FFT pair. Channel identification techniques for other SOPs that exploit 
fi^uency and possibly time correlation in the similar fashion will be obvious to one skilled in 

30 the art. 

It is assumed that a certain subset of available SOP bins are allocated for training. Let 
y be this set of frequency bins used for learning a SIMO channel. To begin, assume that J 
contains v bin indices. Furthermore, let Z, be the v trainmg symbols and X,^ be the 

received data in the training fi^uency-bins from antenna /. Let the quantities , be the 
35 estimated time-domain and frequency-domain channels from the transmit antenna under 
consideration to the receive antenna i. In other words, h, is the v -length impulse response 

from the input under consideration to output /. Likewise, H, is a vector of frequency 
domain values for this channel. With these definitions, it can be shown that 

h,=Y;l(x,,®-'Z,), (81) 

40 and 

H,=Y^^A. (82) 



where. 



v = $[,2,...,i/}. 
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and 



5 This also generalizes to any number of training tones, y > which case the set J includes y 
bin indices. When ;^ ^ v , the frequency domain channel can be determined by, 

= ,{Yj:,Y,,-)r' Y,^,(X,, ®"» Z). (83) 
Note that many of the above calculations can be performed in advance if the training bins are 
predetermined and fixed. Then, the matrix Y^^ .(Y^rYjrjjrV^j. can be computed and 
10 stored. 

Note that there is no requirement that the training symbols always reside in the same 
bins from burst to burst. As long as the transmitter and receiver both know where the 
training symbols are placed in any given burst, the training bins may be varied from one burst 
to the next. This may be usefijl to characterize the nature of colored (across SOP bms) noise 
1 5 and/or interference are present. 

A highly advantageous simplification of (83) can be done when v training sjnnbols are 
placed in bins tfiat are evenly spaced throughout the burst. In other words, 
*/=K^>^>"M^^T^}- In this case, Y;l isequaltothe k -pomt BFFT matrix so that 

S equation (81) represents the execution of an v -point IFFT. One may then obtain in 

20 equation (82) by performing an N-point FFT on a vector consisting of h, padded with A'^ - 

zeros. This approach to identifying H,. is only of computation order {N + v)log2 v . 



Identification over multiple bursts 

Identification accuracy can be improved by either increasing the number of training 
25 symbols within each burst or by averaging over multiple bursts if the channel is correlated 
from one burst to another. Some degree of time domain correlation exists in the channel 
because the Doppler frequency shifts caused by moving objects in the physical environment 
are band-limited. This time correlation can be exploited by recursively filtering the estimated 
channel from the present burst vdth channel estimates from previous bursts. A general 
30 filtering approach is represented by 

h{k + 1) = F(Jfc)h(t) + GQc)h{k) 

where h is the smoothed channel estimate of h over bursts k. The particular recursive filter 
weights F(k) and G(k) can be derived in a number of feshions. Two exemplary filtering 
methods are given in Ae following. The first approach determines a time-invariant FIR filter 
35 for each element of h based on a MMSE cost fimction. The second design is time-varying 
Kabnan filter. 

A particularly simple, yet effective, filter design technique is the determination of a 
time-invariant FIR filter, w, that minimizes the MMSE between the true channel impulse 
response and the filtered estimate. This design approach is referred to as Wiener filtering. In 
40 this embodiment, independent fading is assumed on each element of the channel impulse 
response. Therefore, each element of h can be considered independently. An FIR filter 
produces a filtered estimate by forming a weighted sum of the previous /H-1 estimates for that 
particular impulse response element. 
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10 



15 



20 



25 



Vi = l,2,...,K. 



hXk-p) 

Using V such identical filters for each element of the impulse response, then the 

filtered estimate is given by Yi{k) = (A:) • • • {k)^ , The Wiener filter solution for w 
satisfies the following equation. 



min ^sv"h^ (k) - h, {kf j = min £ 



r 








2\ 




















J 



The solution for the above optimization problem is given by, 

w = [E(h,ik)hXkffE[Mk)hM= • 
If each delay in the channel impulse response undergoes Raleigh &ding is assumed then, 



J,{ap(N+v)T) 



and 



^o(o) 



where T is the sampling rate, q} is the maximum Doppler fi-equency, and Jq is the zeroth- 
order Bessel fimction. The quantities al and are the average channel power and the 
channel estimation noise power, respecttvely. 

This filtering approach has many advantages. First, it is computationally simple. Each 
coefficient of the channel impiilse response is filtered independently with a constant, 
precomputed FIR weighting. Second, the underlying time-correlation in the multipath fading 
channel is efficiently e3q)loited. Third, the exact values used for the filter are optimal in a 
MMSE sense. 

A more generalized time-varying filtering approach is now developed based on the 
Kalman filtering equations. A general model for the time-correlated nature of the channel 
impulse response is given by the following set of equations, 

f(A + l) = Af(jfc)+q(Jfc) 
h(A) = CfW-fr(*) 

where the q and r represent noises with covariances Q and R, respectively. The matrices 
A,C,Q»R are used to define the particular model for the correlation of the impulse response 
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over bursts. Note that the vector h can also include the impluse response coefiScients for 
more than one receive antenna. In this case, the above model can include both time 
correlation and correlation across space. 

In a multi-access scheme, successive channel identifications may occur at an irregular 
5 rate. In this case, this Kahnan filter approach is particularly usefiil since the filtering can be 
done with measurement updates and time updates, 

f (it + 1) = A(I - L(jt)C)f (*) + AL(k)hik) 

Uk) = p(jfc)c^ (cp(jfc)c'' + rY 

P(Jfc) = A(l-L(>t)C)P(Jt)A^ +Q 
h(Jt) = Cf(ft) 

1 0 where L(A) = 0 when the receiver is not receiving data in the present burst. 



Interference Subspace Identification 

For many of the spatial processing techniques embodied in this invention, the 
operation of the TSP and RSP can depend, in part, on the level of interference present in the 
1 5 wireless environment within which the invention operates. More specifically, it may be 

preferable to reduce the amount of mterference contributed to other receivers by a judicious 
choice of the TSWs. It may also be preferable to improve the signal quality at the receiver by 
using RSWs that reject intenrference. In these cases, some quantitative measure of the 
mterference across space and fi^quency is needed. 
20 One preferable measure of the interference present is the so-called interference spatial 

covariance matrix, which describes interference correlation across space for each frequency 
bin, 

R;(«)=4l(«)I(«f], (1) 

where Xj (n) represents an -length received vector of signals fi-om the interfering 

25 transmitter(s). To be more precise, R;(w) describes the interference and noise correlation 
across space for each firequency bin. Since we assume that the noise at the output of each 
receiver antenna path is additive thermal noise, and therefore that the additive noise is 
uncorrected between any two antenna outputs, the noise contribution to R;(«) is non-zero 
only on the matrix diagonal. In environments dominated by interference, i.e. the interference 

30 power at the receiver is much stronger than the additive receiver noise, the noise contribution 
to Rj («) can be neglected. The interference covariance matrix contains information about 
the average spatial behavior of the interference. The eigenvectors of this matrix define the 
average spatial directions (in -space) occupied by the interference. The eigenvalues of 
the matrix indicate the average power occupied by the interference in each the eigendirection. 

35 The eigendirections that are associated with large eigenvalues indicate spatial directions that 
receive a large amount of average interference power. The eigendirections associated with 
small eigenvalues indicate spatial directions that are preferable in that they receive less 
average interference power. 

Identifying the receive covariance matrix, R/(w), is required for finding preferable 

40 RSPs. An analogous transmit covariance matrix is required for finding preferable TSPs. 

Notice that we've defined R/(/j) in terms of received signal samples in Equation (1). Since 
the received signal samples are not usually available at the transmitter, it is preferable to 
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derive the transmit covariance matrix from the receive covariance matrix. In time division 
duplex (TDD) systems, the receive and transmit covariance matrices are substantially equal 
when the time between reception and transmission is short relative to the rate of time 
variation in the channel. In frequency division duplex (FDD) systems, the transmit and 
receive channel values are generally not correlated with one another at any given instant in 
time. However, the transmit and receive covariance matrices are substantially equal in FDD 
systems when sufficient time averaging is used in the calculation of . There are many 
techniques for determining the interference covariance matrix, two of which are discussed 
below. 

One interference characterization approach simply averages the received antenna 
signals during time periods in which the desired transceiver is not transmitting information. 
Since there is no desired signal arriving at the receiver, the interference (and noise) covariance 
is precisely equal to the measured sample covariance matrix, 

R,(«)=R,(n) = ^2^(«'>M"'^')'' • 

In TDD systems, one can make use of "dead-time" to collect samples from the receiver 
during which time no energy from the transmitting end arrives at the receiver. The "dead- 
time" is approximately equal to the round trip propagation delay between the to ends of the 
wireless communications link, and occurs when a transceiver switches from transmission 
mode to reception mode. In the above equation, A, and are the burst indexes 
corresponding to the first and last bursts received during the dead time. Thus, the 
interference covariance can be estimated with no increase in overhead. 

The interference covariance matrices can also be determined while the desired signal is 
being transmitted to the receiver. One approach involves first determining the interference 
signal and subsequently finding the interference signal covariance. The estimated received 
interference is formed by subtracting the estimated desired signal from the total received 
signal, 

i(/i, A) = x(n, k) - H(w,)t)z(w, A) . 
Therefore, once the channel is identified and the information symbols determined, the 
remaining signal is considered to be interference. The interference covariance matrix for bin 
n, averaged over K bursts is given by, 

It is understood that when estimating the covariance matrix, it may be desirable to 
filter the covariance matrix estimates. It may also be advantageous in certain embodiments to 
determine an average interference covariance matrix across SOP bins. For example, within a 
multiple access system bursts may only be received occasionally, making it difficult to acquu'e 
a sufficient number of bursts with which to form an accurate covariance matrix for each bin. 
So instead of averaging over time (a series of received bursts), a covariance matrix is formed 
by averaging over the SOP bins of a single burst, 

^j{k)^jjj;i{n^)i(n,kY , 

It may also be preferable to estimate the interference covariance matrices in an 
alternate frequency band. This can be done using the "dead-time" approach given above. 
This may be advantageous when the transceiver has the capability of choosing alternate 
frequency bands for communicating. Estimates of interference ui alternate bands provides the 
foundation for an adaptive frequency hopped scheme. 
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It is understood that the examples and embodiments described herein are for 
illustrative purposes only and that various modifications or changes in light thereof will be 
suggested to persons skilled in the art and are to be included within the spirit and purview of 
this application and scope of the appended claims and then* fiill scope of equivalents. For 
example, much of the above discussion concerns signal processing in the context of a wireless 
communication system where multiple inputs or multiple outputs are accessed by multiple 
transmitter antenna elements or multiple receiver antenna elements. However, the present 
invention is also usefiil in the context of wireline channels accessible via multiple inputs or 
multiple outputs. 
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1 WHATTS CT.ATMED fS; 

2 1. In a digital communication system, a method for communicating 

3 comprising the steps of: 

4 transmitting signals from one or more transmitter antenna elements; 

5 receiving said signals from via a plurality of receiver antenna elements; 

6 wherein separation of radiation patterns among either said transmitter antenna 

7 elements or said receiver antenna elements is insufficient to establish completely 

8 isolated spatial directions for communication; and wherein 

9 at least one of said transmitting and receiving steps comprises processing said 

10 signals to increase isolation between spatial directions employed for communication at a 

1 1 common frequency. 
12 

13 2. The method of claim 1 wherein a channel coupling said plurality of 

14 transmitter antenna elements and receiver antenna elements at said common frequency is 

15^, characterized by a spatial channel matrix having a rank greater than one. 

r " — 

-16 

17 3. In a digital communication system, a method for communicating 

18 comprising the steps of; 

19 transmitting signals from one or more transmitter antenna elements; 

20 receiving said signals via a plurality of receiver antenna elements; 

21 wherein separation of radiation patterns among either said transmitter antenna 

22 elements or said receiver antenna elements is insufficient to establish completely 

23 isolated spatial directions for communication; and wherein 

24 at least one of said transmitting and receiving steps comprises processing said 

25 signals to increase isolation between subchannels, each subchannel associated with a spatial 

26 direction and a bin of a substantially orthogonalizing procedure. 
27 

28 4. The method of claim 3 wherein said substantially orthogonalizing 

29 procedure belongs to a group including: an inverse Fast Fourier Transform, a Fast Fourier 

30 Transform, a Hilbert transform, a wavelet transform, and processing through a set of bandpass 

3 1 filter/frequency upconverter pairs operating at spaced apart frequencies, . 
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33 5. In a digital communication system, a method for preparing a sequence of 

34 symbols for transmission via a plurality of inputs of a channel: j 

35 a) inputting said symbols of said sequence into a plurality of inputs 

36 corresponding to a plurality of subchannels of said channel, each subchannel corresponding to 

37 an input bin of a transmitter substantially orthogonalizing procedure and a spatial direction; 

38 b) for each input bin, spatially processing symbols inputted to said subchannels 

39 corresponding to said input bin, to develop a spatially processed symbol to assign to each 

40 combination of channel input and input bin of said transmitter substantially orthogonalizing 

41 procedure; and 

42 c) applying, independently for each said channel input, said transmitter 

43 substantially orthogonalizing procedure to said spatially processed symbols assigned to each 

44 said channel input. 
45 

46_ 6. The method of claim 5 wherein said b) step has the effect of making 

47 spatial directions of said subchannels into a set of orthogonal spatial dimensions. 

48 

49 7. The method of claim 5 wherein said transmitter substantially 

50 orthogonalizing procedure belongs to one of a group consisting of an inverse Fast Fourier 

51 Transform, a Fast Fourier Transform, a discrete cosine transform, a Hilbert transform, a 

52 wavelet transform, and processing through a plurality of bandpass filter/frequency converter 

53 pairs centered at spaced apart frequencies. 
54 

55 8. The method of claim 5 further comprising the step of, after said c) step, 

56 applying a cyclic prefix processing procedure to a result of said substantially orthogonalizing 

57 procedure independently for each channel input. 
58 

59 9. The method of claim 5 wherein said transmitter substantially 

60 orthogonalizing procedure is optimized to reduce interference to unintended receivers. 
61 




wo 98/09385 PCT/US97/15160 

39 

62 10. The method of claim 5 wherein said b) step comprises, for each 

63 particular input bin, multiplying a vector comprising symbols allocated to subchannels 

64 corresponding to said input bin by a beneficial weighting matrix, elements of a result vector of 



65 said multiplying step corresponding to different channel inputs of said plurality of channel 

66 inputs. 
67 

68 11. The method of claim 10 wherein said beneficial weighting matrix 

69 comprises an input singular matrix of a matrix containing values representing characteristics of 

70 said- channel, said coupling said plurality of channel inputs to one or more channel outputs. 
71 

72 12. The method of claim 10 wherein said beneficial weighting matrix is 

73 obtained from a matrix containing values representing characteristics of a channel coupling 

74 said plurality of channel inputs to one or more channel outputs. 
75 

76 ^ ^ 13. The method of claim 10 wherein said beneficial weighting matrix is 

77 chosen to reduce interference to unintended receivers. 
78 

79 .14. The method of claim 13 wherein said beneficial weighting matrix is 

80 chosen based upon characterization of a desired signal subspace. 
81 

82 15. The method of claim 14 wherein said beneficial weighting matrix is 

83 chosen fiirther based upon characterization of an imdesired signal subspace. 
84 

85 16. The method of claim 15 wherein characterizations of said desired signal 

86 subspace and said undesired signal subspace are averaged over at least one of time and 

87 frequency. 
88 

89 17. The method of claim 10 wherein said b) step comprises performing said 

90 spatial processing step so as to reduce interference radiated to unintended receivers. 

91 18, The method of claim 10 wherein said b) step comprises, for each input 

92 bin, allocating symbols to each combination of channel input and input bin so that there 
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93 is a one-to-one mapping between spatial direction of a particular subchannel to which a 

94 particular symbol has been allocated and channel input to which said particular symbol 

95 is allocated. 
96 

97 19. The method of claim 10 further comprising the step of prior to said b) 

98 step applying a coding procedure to said symbols. 
99 

100 20. The method of claim 19 wherein said coding procedure is applied 

101 independently for each of said subchannels. 
102 

103 21. The method of claim 19 wherein said coding procedure is applied 

104 independently for each group of subchannels corresponding to an input bin of said substantially 

105 orthogonalizing procedure. 
106 

107^ 22. The method of claim 19 wherein said coding procedure is applied 

independently for each group of subchannels corresponding to a particular spatial direction. 

109 

110 23. The method of claim 19 wherein said coding procedure is applied 

111 integrally across all of said subchaimels. 
112 

113 24. The method of claim 19 wherein said coding procedure belongs to a 

1 14 group consisting of: convolutional coding, Reed-Solomon coding, CRC coding, block coding, 

115 trellis coding, turbo coding, and interleaving. 
116 

117 25. The method of claim 19 wherein said coding procedure comprises a 

118 trellis coding procedure. 
119 

120 26. The method of claim 25 wherein a code design of said trellis coding 

121 procedure is based on one of: improved bit error performance in interference channels, a 

122 periodic product distance metric, exhaustive code polynomial search for favorable bit error 

123 rate polynomial searches, combined weighting of product distance and Euclidean distance. 
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124 product distance of multiple Euclidean distances over short code segments or over a multi- 

125 dimensional symbol, and sum of product distances over short code segments. 
126 

127 27. The method of claim 25 wherein a code design of said trellis coding 

128 procedure is optimized for performance in a fading matrix channel. 
129 

130 28. The method of claim 19 wherein said coding procedure comprises a one- 

131 dimensional trellis coding procedure followed by an interleaving procedure with sequential 

132 groups of symbols output by said trellis coding having their internal order maintained by said 

133 interleaving procedure. 
134 

135 29. The method of claim 19 wherein said coding procedure comprises a 

136 multi-dimensional trellis coding procedure followed by an interleaving procedure with groups 

137 of one-dimensional symbols output simultaneously by said multi-dimensional trellis coding 

138 procedure having their internal order maintained by said interleaving procedure. 



140 30. The method of claim 10 wherein bit loading and power are allocated to 

141 each subchannel. 
142 

143 31. The method of claim 10 further comprising the step of retransmitting 

144 symbols by repeating at least one of said a), b), and c) steps upon receipt of a notification that 

145 said symbols to be retransmitted have been incorrecdy received. 
146 

147 32. The method of claim 10 wherein said channel comprises a wireless 

148 channel and said plurality of channel inputs are associated with a corresponding plurality of 

149 transmitter antenna elements 
150 

151 33. The method of claim 32 wherein said plurality of transmitter antenna 

152 elements are co-located. 
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154 34. The method of claim 32 wherein said plurality of transmitters are at 

155 disparate locations. 
156 

157 35. A method of processing a sequence of symbols received via a plurality of 

158 outputs of a channel, said method comprising the steps of: 

159 a) applying a receiver substantially orthogonalizing procedure to said sequence 

160 of symbols, said procedure being applied independently for each of said plurality of channel 

161 outputs, each output symbol of said receiver substantially orthogonalizing procedure 

162 corresponding to a particular output bin and a particular one of said channel outputs; and 

163 b) for each output bin, spatially processing symbols corresponding to said 

164 output bin to develop spatially processed symbols assigned to a plurality of spatial directions, 

165 each combination of spatial direction and output bin specifying one of a plurality of 

166 subchannels. 
167 

168. 36. The method of claim 35 wherein said b) step has the effect of making 

fef^' said plurality of spatial directions into a set of orthogonal spatial dimensions. 

170 

171 37. The method of claim 35 wherein said receiver substantially 

172 orthogonalizing procedure belongs to one of a group consisting of an inverse Fast Fourier 

173 Transform, a Fast Fourier Transform, a discrete cosine transform, a Hilbert transform, a 

174 wavelet transform, and processing through a plurality of bandpass filter/frequency converter 

175 pairs centered at spaced apart frequencies. 
176 

177 38. The method of claim 35 further comprising the step of, prior to said a) 

178 step, applying a cyclic prefix removal procedure to said sequence of symbols independently 

179 for each of said channel outputs. 
180 

181 39. The method of claim 35 wherein said receiver substantially 

182 orthogonalizing procedure is optimized to reduce deleterious effects of interference from 

183 undesired co-channel transmitters. 
184 
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185 40. • The method of claim 35 wherein said b) step comprises, for each 

186 particular output bin, multiplying a vector comprising symbols of said output bin by a* 

187 beneficial weighting matrix, elements of a result vector of said multiplying step corresponding 

188 to different spatial directions. 
189 

190 41. The method of claim 40 wherein said beneficial weighting matrix 

191 comprises an output singular vector of a matrix containing values representing characteristics 

192 of said channel, said channel coupling one or more channel inputs to said plurality of channel 

193 outputs. 
194 

195 42. The method of claim 40 wherein said beneficial weighting matrix is 

196 chosen to minimize deleterious effects of interference from undesired transmitters, 
197 

198 43, The method of claim 42 wherein said beneficial weighting matrix is 

199. .., chosen based upon characterization of a desired signal subspace. 

201 44. The method of claim 43 wherein said beneficial weighting matrix is 

202 chosen further based upon characterization of an undesired signal subspace. 
203 

204 45. The method of claim 44 wherein said characterizations of said desired 

205 signal subspace and said undesired signal subspace are averaged over at least one of time and 

206 frequency. 
207 

208 46, The method of claim 40 wherein said beneficial weighting matrix is 

209 obtained from a matrix containing values representing characteristics of said charmel, said 

210 channel coupling one or more channel inputs and said plurality of channel outputs. 
211 

212 47. The method of claim 46 wherein said beneficial weighting matrix is 

213 obtained by an MMSE procedure. 
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215 48. The method of claim 35 further comprising the step of after said b) step 

216 applying a decoding procedure to said symbols. 
217 

218 49. The method of claim 48 wherein said decoding procedure is applied 

219 independently for each of said plurality of subchannels. 
220 

221 50. The method of claim 48 wherein said decoding procedure is applied 

222 independently for each group of subchannels corresponding to an output bin of said 

223 substantially orthogonalizing procedure. 
224 

225 51. The method of claim 48 wherein said decoding procedure is applied 

226 independently for each group of subchannels corresponding to a spatial direction. 
227 

228 52. The method of claim 48 wherein said decoding procedure is applied 

229^ ^ integrally across all of said plurality of subchannels, 

231 53. The method of claim 48 wherein said decoding procedure belongs to a 

232 group consisting of Reed-Solomon decoding, CRC decoding, block decoding, and de- 

233 interleaving. 
234 

235 54. The method of claim 48 wherein said decoding procedure comprises a 

236 code sequence detection procedure to decode a trellis code, or convolutional code. 
237 

238 55. The method of claim 54 wherein said code sequence detection procedure 

239 employs a metric belonging to a group consisting of: Euclidean metric, weighted Euclidean 

240 metric, and Hamming metric. 
241 

242 56. The method of claim 48 wherein said decoding procedure reduces 

243 deleterious effects of interference from undesired transmitters. 
244 

245 57. The method of claim 35 further comprising the step of: 
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246 sending a retransmission request when received symbols are 

247 determined to include errors. , 
248 

249 58. The method of claim 35 wherein said channel comprises a wireless 

250 channel and said plurality of channel outputs are coupled to a plurality of corresponding 

251 receiver antenna elements. 
252 

253 59. The method of claim 35 wherein said plurality of receiver antenna 

254 elements are co-located. 
255 

256 60. The method of claim 35 wherein said plurality of receiver antenna 

257 elements are at disparate locations. 
258 

259 61. In a digital communication system, a method for preparing a sequence of 

26(^^ symbols for transmission via a plurality of inputs to a chaimel, said method comprising the 
steps of 

262 selecting a weighting vector for optimal transmission; 

263 applying a transmitter substantially orthogonalizing procedure to 

264 said sequence of symbols to develop a time domain symbol sequence; and 

265 multiplying at least one symbol of said time domain symbol 

266 sequence by said weighting vector to develop a result vector, elements of said result vector 

267 corresponding to symbols to be transmitted via individual ones of said plurality of channel 

268 inputs. 
269 

270 62. The method of claim 61 wherein said weighting vector comprises an 

271 element indicating delay to be applied for a particular one of said plurality of channel inputs. 
272 

273 63. The method of claim 61 wherein said weighting vector is optimized to 

274 reduce interference to unintended receivers. 
275 
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276 64. The method of claim 61 wherein said weighting vector is chosen based 

277 upon characterization of a desired signal subspace. 
278 

279 65. The method of claim 64 wherein said weighting vector is chosen further 

280 based upon characterization of an undesired signal subspace. 
281 

282 66. The method of claim 65 wherein said characterizations of said desired 

283 signal subspace and said undesired signal subspace are averaged over at least one of time and 

284 frequency. 
285 

286 67. The method of claim 61 wherein said channel comprises a wireless 

287 channel and said plurality of channel inputs are associated with a plurality of transmitter 

288 antenna elements. 
289 

290 
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68. In a digital communication system, a method for processing a plurality 

293 of symbols received via a plurality of outputs of a channel, said method comprising the steps 

294 of: 

295 selecting a weighting vector for optimal reception; 

296 multiplying an input vector whose elements correspond to 

297 symbols received substantially simultaneously via a selected one of said plurality of channel 

298 outputs by said weighting vector to obtain a time domain symbol corresponding to a particular 

299 input bin of a receiver substantially orthogonaiizing procedure; 

repeating said multiplying step for successive received symbols to 

301 obtain time domain symbols corresponding to successive input bins of said receiver 

302 substantially orthogonaiizing procedure; and 

applying said receiver substantially orthogonaiizing procedure to 
304 said obtained time domain symbols. 
305 

lP%r 69. The method of claim 68 wherein said weighting vector comprises an 

307 element indicating delay to be applied for a particular one of said plurality of channel outputs. 
308 

309 70. The method of claim 68 wherein said weighting vector is optimized to 

310 reduce deleterious effects of interference from unintended transmitters. 
311 

^^2 71. The method of claim 68 wherein said weighting vector is chosen based 

313 upon characterization of a desired signal subspace. 

314 

^ 72. The method of claim 71 wherein said weighting vector is chosen further 

316 based upon characterization of an undesired signal subspace. 

317 

73. The method of claim 72 wherein said characerizations of said desired 

319 signal subspace and said undesired signal subspace are averaged over at least one of frequency 

320 and time. 
321 
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322 74. The method of claim 71 wherein said chamiel comprises a wireless 

323 chamiel and said plurality of channel outputs are associated with a plurality of corresponding 

324 receiver antenna elements. 
325 

326 75. In a digital communication system, a method of preparing symbols for 

327 transmission via a plurality of inputs of a channel, said method comprising the steps of: 

328 directing symbols to input bins of a transmitter substantially 

329 orthogonalizing procedure so that each input bin has an allocated symbol; 

330 for each particular input bin, spatially processing said symbol 

331 allocated to said particular input bin to develop a spatially processed symbol vector, each 

332 element of said spatially processed symbol vector being assigned to one of said chaimel 

333 inputs; 

334 applying said transmitter substantially orthogonalizing procedure 

335 for a particular channel input, inputs to said substantially orthogonalizing procedure being for 
336^ each input bin, a symbol of said processed symbol vector for said input bin corresponding to 

337 said particular channel input; and 

338 repeating said applying step for each of said plurality of channel 

339 inputs. 
340 

341 76. The method of claim 75 further comprising the step of: 

342 applying a cyclic prefix processing procedure to outputs of said 

343 substantially orthogonalizing procedure independently for each particular channel input. 
344 

345 77. The method of claim 75 wherein said transmitter substantially 

346 orthogonalizing procedure is optimized to reduce interference to unintended receivers. 
347 

348 78. The method of claim 75 wherein said processing step comprises: 

349 multiplying said symbol allocated to said particular input bin by a 

350 beneficial weighting vector to obtain said spatially processed symbol vector. 
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352 79. The method of claim 78 wherein said beneficial weighting vector is an 

3 53 input singular vector of a matrix storing values indicative of said channel, said channel 

354 coupling said plurality of channel inputs and one or more channel outputs. 
355 

356 80. The method of claim 78 wherein said beneficial weighting vector is 

357 chosen to select a beneficial spatial direction for transmission. 
358 

359 81. The method of claim 80 wherein said beneficial weighting vector is 

360 chosen to reduce interference to unintended receivers. 
361 

362 82. The method of claim 81 wherein said beneficial weighting vector is 

363 chosen based upon characterization of a desired signal subspace 
364 

365 83. The method of claim 82 wherein said beneficial weighting vector is 

36d^. chosen fiirther based upon characterization of an undesired signal subspace. 

^557 

368 84. The method of claim 83 wherein said characterizations of said desired 

369 signal subspace and said undesired signal subspace are averaged over at least one of time and 

370 fi-eqiiency. 
371 

372 85. The method of claim 75 wherein said channel comprises a wireless 

373 channel and said plurality of channel inputs are associated with a corresponding plurality of 

374 transmitter antenna elements. 
375 

376 86. In a digital communication system, a method for processing symbols 

377 received by a plurality of outputs of a channel comprising the step of: 

378 applying a receiver substantially orthogonaliang procedure to symbols received 

379 via a particular one of said channel outputs; 

380 repeating said applying step for each of said channel outputs to develop a result 

381 vector for each of a plurality of output bins of said receiver substantially orthogonalizing 

382 procedure, said result vector including a result symbol for each of said channel outputs; and 
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383 for each particular output bin of said receiver substantially orthogonalizing 

384 procedure, spatially processing said result vector for said particular output bin to develop a 

385 spatially processed result symbol for said particular output bin. 
386 

387 87. The method of claim 86 further comprising the step of: 

388 prior to said applying step, applying a cyclic prefix removal procedure to 

389 symbols independently for each. of said channel outputs. 
390 

391 88. The method of claim 86 wherein said substantially orthogonalizing 

392 procedure is optimized to reduce deleterious effects of interference from unintended 

393 transmitters. 
394 

395 89. The method of claim 86 wherein said spatially processing step comprises 

396 multiplying a beneficial weighting vector by said result vector to obtain said spatially 
397^^ processed result symbol. 

398 

399 90. The method of claim 88 wherein said beneficial weighting vector is an 

400 input singular vector of a matrix storing values indicative of characteristics of said channel, 

401 said channel coupling one or more chanel inputs and said plurality of channel outputs. 
402 

403 91, The method of claim 88 wherein said beneficial weighting vector is 

404 chosen to select a particular spatial direction for reception. 
405 

406 92. The method of claim 91 wherein said beneficial weighting vector is 

407 chosen to minimize deleterious effects of interference fi-om unintended transnutters. 
408 

409 93. The method of claim 91 wherein said beneficial weighting vector is 

410 chosen based upon characterization of a desired signal subspace. 
411 

412 94. The method of claim 93 wherein said beneficial weighting vector is 

413 chosen based upon characterization of an undesired signal subspace. 




wo 98/09385 PCT/US97/15160 

51 



414 

415 95. The method of claim 94 wherein said characterizations of said desired 

416 signal subspace and said undesired signal subspace are averaged over at least one of time and 

417 frequency. 
418 

419 96. The method of claim 86 wherein said channel comprises a wireless 

420 channel and said plurality of channel outputs are associated with a corresponding plurality of 

421 channel outputs. 
422 

423 97. In a digital communication system including a communication channel 

424 having one or more inputs and at least one or more outputs, a method for determining 

425 characteristics of said channel based on signals received by said one or more outputs, 

426 comprising the steps of: 

427 a) receiving via said one or more channel outputs, at least v training symbols 
428^.. transmitted via a particular spatial direction of said channel, v being an extent in symbol 
^9 periods of a duration of significant terms of an impulse response of a channel; and 

430 b) applying a substantially orthogonalizing procedure to said received at least 

43 1 vtraining symbols to obtain a time domain response for said spadal direction; and 

432 c) applying an inverse of said substantially orthogonaliang procedure to a zero- 

433 padded version of said time domain response to obtain a frequency response for said particular 

434 spatial direction. 
435 

436 98. The method of claim 97 wherein said substantially orthogonalizing 

437 procedure comprises an inverse Fast Fourier Transform and said inverse of said substantially 

438 orthogonalizing procedure comprises a Fast Fourier Transform. 
439 

440 99. The method of claim 98 wherein said a) step comprises receiving exactiy 

441 V training symbols. 
442 

443 100. The method of claim 97 further comprising the step of repeating said a), 

444 b), c), and d) steps for a plurality of spatial directions. 
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445 

446 101. The method of claim 99 wherein each of said plurality of spatial 

447 directions corresponds to transmission through one of said plurality of channel inputs 

448 exclusively. 
449 

450 102. The method of claim 98 wherein said v training symbols belong to a 

451 burst of N symbols and said characteristics are determined for said burst. 
452 

453 103. The method of claim 102 further comprising the steps of repeating said 

454 a), b), c), and d) steps for successive bursts. 
455 

456 104. The method of claim 103 further comprising the step of after, said b) 

457 step, smoothing said time-domain response over successive bursts. 
458 

459 105. The method of claim 104 wherein said smoothing step comprises Kalman 

^^^^ filtering. 

461 

462 106. The method of claim 104 wherein said smoothing step comprises Wiener 

463 filtering. 
464 

465 107. The method of claim 97 wherein said communication channel comprises 

466 known and unknown components, wherein said effects of said known components are removed 

467 by deconvolution, and characteristics of said unknown components are determined by said a), 

468 b), c), and d) steps, thereby reducing . 
469 

470 108. In a digital communication system including a communication channel 

471 having one or more inputs and one or more outputs, a method for determining characteristics 

472 of said channel based on signals received via one or more channel outputs, comprising the 

473 steps of: 

474 receiving training symbols via said channel outputs; and 
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475 computing characteristics of said channel based on said received 

476 training, symbols and assumptions that an impulse response of said channel is substantially 

477 time-limited and that variation of said impulse response over time is continuous. 
478 

479 109. In a digital communication system, a method for communicating over a 

480 channel having at least one input and at least one output, and having a plurality of either inputs 

481 or outputs, said method comprising the steps of: 

482 dividing said channel into a plurality of subchannels, each 

483 subchannel corresponding to a combination of spatial direction and an input bin of a 

484 substantially orthogonalizing procedure; and 

485 communicating symbols over one or more of said plurality of 

486 subchannels. 
487 

488 110. In a digital communication system, a method for preparing a sequence of 

489^*, symbols for transmission via a plurality of inputs of a channel, comprising the steps of: 

490 a) inputting said symbols of said sequence into a plurality of 

491 input corresponding to a plurality of subchannels of said channel, each subchannel 

492 corresponding to an input bin of a transmitter substantially orthogonalizing procedure and a 

493 channel input; and 

494 b) applying, independentiy for each said chaimel input, said 

495 transmitter substantially orthogonalizing procedure to said symbols assigned to each said 

496 channel input. 
497 

498 111. A method of processing a sequence of symbols received via a plurality of 

499 outputs of a channel, said method comprising the steps of: 

500 a) applying a substantially orthogonalizing procedure to said 

501 sequence of symbols, said procedure being applied independentiy for each of said plurality of 

502 channel outputs, each output symbol of said substantially orthogonalizing procedure 

503 corresponding to a subchannel identified by a combination of a particular output bin and a 

504 particular one of said channel outputs; and 

505 b) processing symbols in said subchannels. 
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112. In a digital communication system, apparatus for communicating 



compnsmg: 



antenna elements; 
receiver antenna elements; 



a transmitter that transmits signals from one or more transmitter 
a receiver that receives said signals from via a plurality of 



wherein separation of radiation patterns among either said 
transmitter antenna elements or said receiver antenna elements is insufficient to 
establish completely isolated spatial directions for communication; and wherein 

at least one of said transmitter and said receiver comprises a 
processor that processes said signals to increase isolation between spatial directions employed 
for communication at a common frequency. 

113. The apparatus of claim 1 12 wherein a channel coupling said plurality of 
transmitter antenna elements and receiver antenna elements at said common frequency is 
characterized by a spatial channel matrix having a rank greater than one. 



114. In a digital communication system, apparatus for communicating 



compnsmg: 



antenna elements; 
antenna elements; 



a transmitter transmitting signals from one or more transmitter 
a receiver receiving said signals via a plurality of receiver 



wherein separation of radiation patterns among either said 
transmitter antenna elements or said receiver antenna elements is insufficient to 
establish completely isolated spatial directions for communication; and wherein 

at least one of said transmitter and said receiver comprises a 
processor that processes said signals to increase isolation between subchannels, each 
subchannel associated with a spatial direction and a bin of a substantially orthogonalizing 
procedure. 
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537 

538 115. The apparatus of claim 1 14 wherein said substantially orthogonalizing 

539 procedure belongs to a group including: an inverse Fast Fourier Transform, a Fast Fourier 

540 Transform, a Hilbert transform, a wavelet transform, and processing through a set of bandpass 

54 1 filter/frequency upconverter pairs operating at spaced apart frequencies. . 
542 

543 116. In a digital communication system, apparatus for preparing a sequence of 

544 symbols for transmission via a plurality of inputs of a channel: 

545 a plurality of parallel subchaimel inputs receiving said symbols, said parallel 

546 subchaimel inputs corresponding to a plurality of subchannels, each subchannel corresponding 

547 to an input bin of a transmitter substantially orthogonalizing procedure and a spatial direction; 

548 a spatial processor that, for each input bin, spatially processor symbols received 

549 by said subchannel inputs corresponding to said input bin, to develop a spatially processed 

550 symbol to assign to each combination of channel input and input bin of said transmitter 
551^^, substantially orthogonalizing procedure; and 

532 a substantially orthogonal procedure processor system that applies, 

553 independently for each said channel input, said transmitter substantially orthogonalizing 

554 procedure to said spatially processed symbols assigned to each said channel input. 
555 

556 117. The apparatus of claim 1 16 wherein said spatial processor has the effect 

557 of making spatial directions of said subchannels into a set of orthogonal spatial dimensions. 
558 

559 1 18. The apparatus of claim 116 wherein said transmitter substantially 

560 orthogonalizing procedure belongs to one of a group consisting of an inverse Fast Fourier 

561 Transform, a Fast Fourier Transform, a discrete cosine transform, a Hilbert transform, a 

562 wavelet transform, and processing through a plurality of bandpass filter/frequency converter 

563 pairs centered at spaced apart frequencies. 
564 

565 119. The apparatus of claim 116 further comprising: a cyclic prefix processor 

566 that applies a cyclic prefix processing procedure to a result of said substantially 

567 orthogonalizing procedure independently for each channel input. 
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568 

569 120. The apparatus of claim 1 16 wherein said transmitter substantially 

570 orthogonalizing procedure is optimized to reduce interference to unintended receivers. 
571 

572 121. The apparatus of claim 1 16 wherein said spatial processor comprises, for 

573 each particular input bin, a weight multiplier that multiplies a vector comprising symbols 

574 allocated to subchannels corresponding to said input bin by a beneficial weighting matrix, 

575 elements of a result vector of said weight multiplier corresponding to different channel inputs 

576 of said plurality of channel inputs. 
577 

578 122. The apparatus of claim 121 wherein said beneficial weighting matrix 

579 comprises an input singular matrix of a matrix containing values representing characteristics of 

580 said channel, said channel coupling said plurality of channel inputs to one or more channel 

581 outputs. 

582 ^ 

123. The apparatus of claim 121 wherein said beneficial weighting matrix is 

584 obtained from a matrix containing values representing characteristics of a channel coupling 

585 said plurality of channel inputs to one or more channel outputs. 
586 

587 124. The apparatus of claim 121 wherein said beneficial weighting matrix is 

588 chosen to reduce interference to unintended receivers. 
589 

590 125. The apparatus of claim 124 wherein said beneficial weighting matrix is 

591 chosen based upon characterization of a desired signal subspace. 
592 

593 126. The apparatus of claim 125 wherein said beneficial weighting matrix is 

594 chosen fijrther based upon characterization of an undesired signal subspace. 
595 

596 127. The apparatus of claim 126 wherein characterizations of said desired 

597 signal subspace and said undesired signal subspace are averaged over at least one of time and 

598 frequency. 
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599 

600 128. The apparatus of claim 1 16 wherein said spatial processor operates so as 

601 to reduce interference radiated to unintended receivers. 
602 

603 129. The apparatus of claim 1 16 wherein said spatial processor, allocates 

604 symbols to each combination of channel input and input bin so that there is a one-to-one 

605 mapping between spatial direction of a particular subchannel to which a particular symbol has 

606 been allocated and channel input to which said particular symbol is allocated. 
607 

608 130. The apparatus of claim 116 further comprising a coder that applies a 

609 coding procedure to said symbols prior to processing by said spatial processor. 
610 

611 131. The apparatus of claim 130 wherein said coding procedure is applied 

612 independently for each of said subchannels. 

613^^ 

132. The apparatus of claim 130 wherein said coding procedure is applied 

615 independently for each group of subchannels corresponding to an input bin of said substantially 

616 orthogonalizing procedure. 
617 

618 133. The apparatus of claim 130 wherein said coding procedure is applied 

619 independentiy for each group of subchannels corresponding to a particular spatial direction. 
620 

621 134. The apparatus of claim 130 wherein said coding procedure is applied 

622 integrally across ail of said subchannels. 
623 

624 135, The apparatus of claim 130 wherein said coding procedure belongs to a 

625 group consisting of: convolutional coding, Reed-Solomon coding, CRG coding, block coding, 

626 trellis coding, turbo coding, and interleaving. 
627 

628 136. The apparatus of claim 130 wherein said coding procedure comprises a 

629 trellis coding procedure. 
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630 

631 137. The apparatus of claim 136 wherein a code design of said trellis coding 

632 procedure is based on one of: improved bit error performance in interference channels, a 

633 periodic product distance metric, exhaustive code polynomial search for favorable bit error 

634 rate polynomial searches, combined weighting of product distance and Euclidean distance, 

635 product distance of multiple Euclidean distances over short code segments or over a multi- 

636 dimensional symbol, and sum of product distances over short code segments. 
637 

638 138, The apparatus of claim 136 wherein a code design of said trellis coding 

639 procedure is optimized for performance in a fading matrix channel. 
640 

641 139. The apparatus of claim 130 wherein said coding procedure comprises a 

642 one-dimensional trellis coding procedure followed by an interleaving procedure with sequential 

643 groups of symbols output by said trellis coding having their internal order maintained by said 
644. interleaving procedure. 

646 140. The apparatus of claim 130 wherein said coding procedure comprises a 

647 multi-dimensional trellis coding procedure followed by an interleaving procedure with groups 

648 of one-dimensional symbols output simultaneously by said multi-dimensional trellis coding 

649 procedure having their internal order mmntained by said interleaving procedure. 
650 

651 141. The apparatus of claim 130 wherein bit loading and power are allocated 

652 to each subchannel. 
653 

654 142. The apparatus of claim 1 16 further comprising an ARQ system that 

655 retransmits symbols via at least one of said spatial processor, and said substantially 

656 orthogonalizing procedure processor upon receipt of a notification that said symbols to be 

657 retransmitted have been incorrectly received. 
658 
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659 143. The apparatus of claim 1 16 wherein said channel comprises a wireless 

660 channel and said plurality of channel inputs are associated with a corresponding plurality of 

661 transmitter antenna elements 
662 

663 144, The apparatus of claim 142 wherein said plurality of transmitter antenna 

664 elements are co-located. 
665 

666 145. The apparatus of claim 144 wherein said plurality of transmitters are at 

667 disparate locations. 
668 

669 146. Apparatus of processing a sequence of symbols received via a plurality 

670 of outputs of a channel, said apparatus comprising: 

671 a substantially orthogonalizing procedure processor system that applies a 

672 receiver substantially orthogonalizing procedure to said sequence of symbols, said procedure 
67^^L being applied independently for each of said plurality of channel outputs, each ou^ut symbol 



of said substantially orthogonalizing procedure corresponding to a particular output bin and a 

675 particular one of said channel outputs; and 

676 a spatial processor that, for each output bin, spatially processes symbols 

677 corresponding to said output bin to develop spatially processed symbols assigned to a plurality 

678 of spatial directions, each combination of spatial direction and output bin specifying one of a 

679 plurality of subchannels. 
680 

681 147, The apparatus of claim 146 wherein said spatial processor operates to 

682 make said plurality of spatial directions into a set of orthogonal spatial dimensions. 
683 

684 148. The apparatus of claim 146 wherein said receiver substantially 

685 orthogonalizing procedure belongs to one of a group consisting of an inverse Fast Fourier 

686 Transform, a Fast Fourier Transform, a discrete cosine transform, a Hilbert transform, a 

687 wavelet transform, and processing through a plurality of bandpass filter/frequency converter 

688 pairs centered at spaced apart frequencies. 
689 
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690 149. The apparatus of claim 146 further comprising: a cyclic prefix processor 

691 that applies a cyclic prefix removal procedure to said sequence of symbols independently for 

692 each of said channel outputs. 
693 

694 150. The apparatus of claim 146 wherein said receiver substantially 

695 orthogonalizing procedure is optimized to reduce deleterious effects of interference from 

696 undesired co-channel transmitters. 
697 

698 151. The apparatus of claim 146 wherein said spatial processor comprises, for 

699 each particular output bin, a weight multiplier that multiplies a vector comprising symbols of 

700 said output bin by a beneficial weighting matrix, elements of a result vector of said multiplier 

701 corresponding to different spatial directions. 
702 

703 152. The apparatus of claim 151 wherein said beneficial weighting matrix 

K comprises an output singular vector of a matrix containing values representing characteristics 

/05 of said chaimel, said channel coupling one or more channel inputs to said plurality of channel 

706 outputs. 
707 

708 153. The apparatus of claim 151 wherein said beneficial weighting matrix is 

709 chosen to minimize deleterious effects of interference from undesired transmitters. 
710 

711 154. The apparatus of claim 151 wherein said beneficial weighting matrix is 

712 chosen based upon characterization of a desired signal subspace. 
713 

714 155. The apparatus of claim 154 wherein said beneficial weighting matrix is 

715 chosen fiuther based upon characterization of an undesired signal subspace. 
716 

717 156. The apparatus of claim 155 wherein said characterizations of said desired 

718 signal subspace and said undesired signal subspace are averaged over at least one of time and 

719 frequency. 
720 
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721 157. The apparatus of claim 151 wherein said beneficial weighting matrix is 

722 obtained from a matrix containing values representing characteristics of said channel, said 

723 channel coupling one or more channel inputs and said plurality of channel outputs. 
724 

725 158. The apparatus of claim 157 wherein said beneficial weighting matrix is 

726 obtained by an MMSE procedure. 
727 

728 159. The apparatus of claim 146 further comprising: a decoder that applies a 

729 decoding procedure to said spatially processed symbols. 
730 

731 160. The apparatus of claim 159 wherein said decoding procedure is applied 

732 independently for each of said plurality of subchannels. 
733 

734 161 xhe apparatus of claim 159 wherein said decoding procedure is applied 

735,. independenfly for each group of subchannels corresponding to an output bin of said 

736 substantially orthogonalizing procedure. 
737 

738 162. The apparatus of claim 159 wherein said decoding procedure is applied 

739 independently for each group of subchannels corresponding to a spatial direction. 
740 

741 163. The apparatus of claim 159 wherein said decoding procedure is applied 

742 integrally across all of said plurality of subchannels. 
743 

744 164. The apparatus of claim 159 wherein said decoding procedure belongs to 

745 a group consisting of: Reed-Solomon decoding, CRC decoding, block decoding, and de- 

746 interleaving. 
747 

748 165. The apparatus of claim 159 wherein said decoding procedure comprises a 

749 code sequence detection procedure to decode a trellis code, or convolutional code. 
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751 166. The apparatus of claim 165 wherein said code sequence detection 

752 procedure employs a metric belonging to a group consisting of: Euclidean metric, weighted 

753 Euclidean metric, and Hamming metric. 
754 

755 167. The apparatus of claim 159 wherein said decoding procedure reduces 

756 deleterious effects of interference from undesired transmitters. 
757 

758 168. The apparatus of claim 146 further comprising; 

759 a system' that sends a retransmission request when received symbols are 

760 determined to include errors. 
761 

762 169. The apparatus of claim 170 wherein said channel comprises a wireless 

763 channel and said plurality of channel outputs are coupled to a plurality of corresponding 

764 receiver antenna elements. 

76^.^ 

766 171. The apparatus of claim 170 wherein said plurality of receiver antenna 

767 elements are co-located. 
768 

769 172. The apparatus of claim 170 wherein said plurality of receiver antenna 

770 elements are at disparate locations. 
771 

772 173. In a digital conmiunication system, apparatus for preparing a sequence of 

773 symbols for transmission via a plurality of inputs to a channel, said apparatus comprising: 

774 a substantially orthogonal procedure processor that applies a transmitter 

775 substantially orthogonalizing procedure to said sequence of symbols to develop a time domain 

776 symbol sequence; and 

777 a weight multiplier that multiplies at least one symbol of said time domain 

778 symbol sequence by a weighting vector selected for improved communication to develop a 

779 result vector, elements of said result vector corresponding to symbols to be transmitted via 

780 individual ones of said plurality of channel inputs. 
781 
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782 174. The apparatus of claim 173 wherein said weighting vector comprises an 

783 element indicating delay to be applied for a particular one of said plurality of channel inputs. 
784 

785 175. The apparatus of claim 174 wherein said weighting vector is optimized 

786 to reduce interference to unintended receivers. 
787 

788 176. The apparatus of claim 173 wherein said weighting vector is chosen 

789 based upon characterization of a desired signal subspace. 
790 

791 177. The apparatus of claim 176 wherein said weighting vector is chosen 

792 further based upon characterization of an undesired signal subspace. 
793 

794 178. The apparatus of claim 177 wherein said characterizations of said desired 

795 signal subspace and said undesired signal subspace are averaged over at least one of time and 



-79&*^. frequency. 



797 

798 179. The apparatus of claim 173 wherein said channel comprises a wireless 

799 channel and said plurality of channel inputs are associated with a 

800 plurality of transmitter antenna elements. 

801 180. 

802 180. In a digital communication system, apparatus for processing a plurality 

803 of symbols received via a plurality of outputs of a channel, said apparatus comprising: 

804 a weight multiplier that performs a multiplication of an input vector whose 



805 elements correspond to symbols received substantially simultaneously via a selected one of said 

806 plurality of channel outputs by a weighting vector to obtain a time domain symbol 

807 corresponding to a particular input bin of a receiver substantially orthogonalizing procedure 

808 and that repeats said multiplication for successive received symbols to obtain time domain 

809 symbols corresponding to successive input bins of said receiver substantially orthogonalizing 

810 procedure; and 

811 a substantial orthogonalizing procediu-e processor that applies said substantially 

812 orthogonalizing procedure processor to said obtained time domain symbols. 
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813 

814 181 . The apparatus of claim 180 wherein said weighting vector comprises an 

815 element indicating delay to be applied for a particular one of smd plurality of channel outputs. 
816 

817 182. The apparatus of claim 180 wherein said weighting vector is optimized 

818 to reduce deleterious effects of interference from unintended transmitters. 
819 

820 183. The apparatus of claim 180 wherein said weighting vector is chosen 

821 based upon characterization of a desired signal subspace. 

822 . 

823 184. The apparatus of claim 183 wherein said weighting vector is chosen 

824 further based upon characterization of an undesired signal subspace. 
825 

826 185. The apparatus of claim 184 wherein said characterizations of said desired 

signal subspace and said undesired signal subspace are averaged over at least one of frequency 

$28 and time. 
829 

830 186. The apparatus of claim 180 wherein said channel comprises a wireless 

831 channel and said plurality of channel outputs are associated with a plurality of corresponding 

832 receiver antenna elements. 
833 

834 187. In a digital communication system, apparatus for preparing s)mibols for 

835 transmission via a plurality of inputs of a chaimel, said apparatus comprising; 

836 a plurality of symbol inputs, each of said symbol inputs receiving a symbol 

837 intended for a particular input bin of a transmitter substantially orthogonalizing procedure so 

838 that each of a plurality of input bins of said transmitter substantially orthongonalizing 

839 procedure has an allocated symbol; 

840 a spatial processor that, for each particular input bin of said plurality of input 

841 bins, spatially processes said symbol allocated to said particular input bin to develop a spatially 

842 processed symbol vector, each element of said spatially processed symbol vector being 

843 assigned to one of said channel inputs; and 
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844 a substantially orthogonalizing procedure processor that applies said 

845 substantially orthogonalizing procedure for a particular channel input, inputs to said 

846 substantially orthogonalizing procedure being for each input bin, a symbol of said processed 

847 symbol vector for said input bin corresponding to said particular channel input, and that 

848 applies said sustantially orthogonalizing procedure for each of said plurality of chaimel inputs. 
849 

850 188. The apparatus of claim 187 further comprising: 

851 a cyclic prefix processor that applies a cyclic prefix processing procedure to 

852 outputs of said substantially orthogonaliang procedure independently for each particular 

853 channel input. 
854 

855 189, The apparatus of claim 187 wherein said substantially orthogonalizing 

856 procedure is optimized to reduce interference to unintended receivers. 
857 

858^. 190. The apparatus of claim 187 wherein said spatial processor comprises: 

a weight multiplier that multiplies said symbol allocated to said particular input 

860 bin by a beneficial weighting vector to obtain said spatially processed sjrmbol vector. 
861 

862 191. The apparatus of claim 190 wherein said beneficial weighting vector is 

863 an input singular vector of a matrix storing values indicative of characteristics of said channel, 

864 said channel coupling said plurality of channel inputs and one or more channel outputs. 
865 

866 192. The qsparatus of claim 190 wherein said beneficial weighting vector is 

867 chosen to select a beneficial spatial direction for transmission. 
868 

869 193. The apparatus of claim 191 wherein said beneficial weighting vector is 

870 chosen to reduce interference to unintended receivers. 
871 

872 194. The apparatus of claim 193 wherein said beneficial weighting vector is 

873 chosen based upon characterization of a desired signal subspace 
874 
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875 195. The apparatus of claim 194 wherein said beneficial weighting vector is 

876 chosen further based upon characterization of an undesired signal subspace. 

877 • 

878 196. The apparatus of claim 195 wherein said characterizations of said desired 

879 signal subspace and said undesired signal subspace are averaged over at least one of time and 

880 frequency. 
881 

882 197. The apparatus of claim 187 wherein said channel comprises a wireless 

883 channel and said plurality of channel inputs are associated with a corresponding plurality of 

884 transmitter antenna elements. 
885 

886 198. In a digital communication system, apparatus for processing symbols 

887 received by a plurality of outputs of a channel comprising: 

888 a substantially orthogonalizing procedure processor that applies a receiver 
889-^- substantially orthogonalizing procedure to symbols received via a particular one of said 

o90 channel outputs and that said applies said receiver substantially orthogonalizing procedure for 

891 each of said channel outputs to develop a result vector for each of a plurality of output bins of 

892 said substantially orthogonalizing procedure, said result vector including a result symbol for 

893 each of said channel outputs; and 

894 a spatial processor that, for each particular output bin of said substantially 

895 orthogonalizing procedure, spatially processes said result vector for said particular output bin 

896 to develop a spatially processed result symbol for said particular output bin. 
897 

898 199. The apparatus of claim 198 ftirther comprising: a cyclic prefix removal 

899 processor that applies a cyclic prefix removal procedure to symbols independentiy for each of 

900 said channel outputs. 
901 

902 200. The apparatus of claim 198 wherein said substantially orthogonaIi2ing 

903 procedure is optimized to reduce deleterious effects of interference from unintended 

904 transmitters. 
905 
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906 201, The apparatus of claim 198 wherein said spatially processor comprises a 

907 weight multiplier that multiplies a beneficial weighting vector by said result vector to obtain 

908 said spatially processed result symbol. 
909 

202. The apparatus of claim 201 wherein said beneficial weighting vector is 

911 an input singular vector of a matrix storing values indicative of characteristics of said channel, 

912 said channel coupling one or more chanel inputs and said plurality of channel outputs. 
913 

914 203. The apparatus of claim 201 wherein said beneficial weighting vector is 

915 chosen to select a particular spatial direction for reception. 
916 

917 204, The apparatus of claim 203 wherein said beneficial weighting vector is 

918 chosen to minimize deleterious effects of interference from unintended transmitters. 
919 

920 ^ 205. The apparatus of claim 204 wherein said beneficial weighting vector is 
chosen based upon characterization of a desired signal subspace. 

922 

923 206. The apparatus of claim 205 wherein said beneficial weighting vector is 

924 chosen based upon characterization of an undesired signal subspace. 
925 

926 207. The apparatus of claim 206 wherein said characterizations of said desired 

927 signal subspace and said undesired signal subspace are averaged over at least one of time and 

928 fi-equency. 
929 

930 208. The apparatus of claim 198 wherein said channel comprises a wireless 

931 channel and said plurality of channel outputs are associated with a corresponding plurality of 

932 channel outputs. 
933 

934 209. In a digital communication system including a communication channel 

935 having one or more inputs and at least one or more outputs apparatus for determining 
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936 characteristics of said channel based on signals received by said one or more outputs, 

937 comprising: 

938 a receiver system receiving via said one or more channel outputs, at least 

939 training symbols transmitted via a particular spatial direction of said channel, being an extent 

940 in symbol periods of a duration of significant terms of an impulse response of a channel; 

941 a substantially orthogonalizing procedure processor that applies a substantially 

942 orthogonalizing procedure processor to said received at least training symbols to obtain a 

943 time domain response for said particular spatial direction; and 

944 an inverse substantially orthogonalizing procedure processor that applies an 

945 inverse of said substantially orthogonalizing procedure to a zero-padded version of said time 

946 domain response to obtain a frequency response for said particular spatial direction. 
947 

948 210. The apparatus of claim 209 wherein said substantially orthogonalizing 

949 procedure comprises an inverse Fast Fourier Transform and said inverse of said substantially 
950_ ^ orthogonalizing procedure comprises a Fast Fourier Transform. 

952 21 1. The apparatus of cldxm 209 wherein said receiver system receives exactiy 

953 training symbols. 
954 

955 212. The apparatus of claim 209 wherein said receiver system, said 

956 substantially orthogonalizing procedure processor and said inverse substantially 

957 orthogonalizing procedure process operate repeatedly for a plurality of spatial directions. 
958 

959 213. The apparatus of claim 209 wherein each of said plurality of spatial 

960 directions corresponds to transmission through one of said plurality of channel inputs 

961 exclusively. 
962 

963 214. The apparatus of claim 209 wherein said training symbols belong to a 

964 burst of N symbols and said characteristics are determined for said burst 
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966 215. The apparatus of claim 214 said receiver system, said substantially 

967 orthogonalizing procedure processor and said inverse substantially orthogonaiizing procedure 

968 process operate repeatedly for a plurality of bursts. 
969 

970 216. The apparatus of claim 215 further comprising: 

971 means for smoothing said time-domain response over successive bursts. 
972 

973 217. The apparatus of claim 216 wherein said smoothing means comprises: 

974 means for Kalman filtering said time-domain response over successive bursts. 
975 

976 218. The apparatus of claim 217 wherein said smoothing means comprises 

977 means for Wiener filtering said time-domain response over successive bursts. 
978 

979 219, The apparatus of claim 209 wherein said communication channel 

98^^ comprises known and unknown components, wherein said effects of said known components 

f S i are removed by deconvolution, and characteristics of said unknown components are 

982 determined by said a), b), c), and d) steps, thereby reducing . 
983 

984 220. In a digital communication system including a conmiimication channel 

985 having one or more inputs and one or more outputs, apparatus for determining characteristics 

986 of said channel based on signals received via one or more channel outputs, comprising: 

987 a receiver that receives training symbols via said channel outputs; and 

988 a processor that computes characteristics of said channel based on said received 

989 training symbols and assumptions that an impulse response of said channel is substantially 

990 time-limited and that variation of said impulse response over time is continuous. 
991 

992 221. In a digital communication system, apparatus for communicating over a 

993 channel having at least one input and at least one output, and having a plurality of either inputs 

994 or outputs, said apparatus comprising: 
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means for dividing said channel into a plurality of subchannels, each subchannel 
corresponding to a combination of spatial direction and an input bin of a substantially 
orthogonalizing procedure; and 

means for communicating symbols over one or more of said plurality of 

subchannels. 

222. In a digital communication system, apparatus for preparing a sequence of 
symbols for transmission via a plurality of inputs of a channel, said apparatus comprising: 

a plurality of parallel subchannel inputs that receive said sequence of symbols, 
said subchannel inputs corresponding to a plurality of subchannels, each subchannel 
corresponding to an input bin of a transmitter substantially orthogonalizing procedure and a 
channel input; and 

a substantially orthogonalizing procedure processor that applies, independently 
for each said channel input, said transmitter substantially orthogonalizing procedure to said 
symbols assigned to each said channel input. 

223. Apparatus for processing a sequence of symbols received via a plurality 
of outputs of a channel, said apparatus comprising the steps of: 

a substantially orthogonalizing procedure processor that applies a receiver 
substantially orthogonalizing procedure to said sequence of symbols, said procedure being 
applied independently for each of said plurality of channel outputs, each output symbol of said 
receiver substantially orthogonalizing procedure corresponding to a subchannel identified by a 
combination of a particular output bin and a particular one of said channel outputs; and 

a processor that processes symbols in said subchannels. 
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